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1 1 Abstract 

CD . 

! Combining dispersion and operator product expansion techniques, we derive the conformal partial 

> ■ 

V~l ' wave decomposition of the virtual Compton scattering amplitude in terms of complex conformal 
spin to twist-two accuracy. The perturbation theory predictions for the deeply virtual Comp- 

c3 ; 

ton scattering (DVCS) amplitude are presented in next-to-leading order for both conformal and 
modified minimal subtraction scheme. Within a conformal subtraction scheme, where we exploit 
predictive power of conformal symmetry, the radiative corrections are presented up to next-to- next- 
to-leading order accuracy. Here, because of the trace anomaly, the mixing of conformal moments 
of generalized parton distributions (GPD) at the three-loop level remains unknown. Within a new 
proposed parameterization for GPDs, we then study the convergence of perturbation theory and 
demonstrate that our formalism is suitable for a fitting procedure of DVCS observables. 
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1 Introduction 



The partonic content of the nucleon, other hadrons, and nuclei has been studied within hard 
inclusive processes over almost four decades. Since partons are confined by the strong force, 
their response to, e.g., an electromagnetic probe cannot be directly accessed in experiments; the 
long-range interaction in their environment remains essential. Fortunately, based on factorization 
theorems [1], the short- and the long-distance interaction can often be separated. The former 
depends on the specific process and can be systematically computed using perturbation theory, 
while the latter is encoded in process-independent nonperturbative quantities, e.g., parton den- 
sities. These densities have a semiclassical interpretation within the parton picture, which is, 
however, not independent of the conventions used in the evaluation of short-distance physics. To- 
gether with an increasing amount of experimental data, the perturbative factorization approach 
leads to a deeper and more precise insight into the hadronic world, which is essential even for the 
search of new physics. 

Most importantly, the factorization approach is a tool that relates observables measured in 
various inclusive processes and thus has a predictive power. The improvement of quantitative 
predictions requires, besides precise experimental measurements, also a refinement of the theo- 
retical approach. This includes both the perturbative evaluation of radiative corrections to the 
short-distance physics at higher orders and an understanding of the so-called power suppressed 
corrections, which alter the factorization theorems. In particular, the effort in the perturbative 
sector, which reached the three-loop level, led to a quantitative understanding of the inclusive 
QCD physics at the level of a few percent. In practice, parton densities have for a long time 
been extracted via a theoretically motivated functional ansatz, see, e.g., Ref. [2], depending on 
a number of parameters, that is globally fitted to experimental data [3, 4, 5, 6, 7, 8, 9, 10]. In 
such fits, based on traditional \ 2 minimization, the error estimation using the hypothesis of linear 
error propagation has become standard in the last few years [3, 4, 5, 8]. To overcome drawbacks 
of this traditional method, alternative frameworks of statistical inference [11] and neural network 
parameterization [12, 13] have been proposed. 

In contrast to inclusive processes, the understanding and the theoretical description of exclusive 
ones remain poor. In general, there is a lack of experimental data and so the underlying theoretical 
framework, e.g., the factorization of hadronic amplitudes in terms of distribution amplitudes 
[14, 15, 16, 17, 18], cannot be quantitatively tested. Moreover, the theoretical framework at the 
next-to-leading order (NLO) is developed for only a few processes and the applicability of the 
factorization approach at accessible scales is controversially discussed in the literature [19, 20, 21]. 
The old problem of whether all valence partons take part in the short-distance process or only the 
active one (Feynman mechanism) remains open. A way out of this theoretical dilemma might be 
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provided by the use of light-cone sum rules [22, 23]. 

Some years ago a new nonperturbative generalized distribution amplitude 1 was proposed as a 
means to access the partonic content of hadrons [24, 25] . In particular, it was theoretically studied 
in connection with deeply virtual Compton scattering (DVCS), where the partonic content of a 
hadron is probed by two photons. To leading order, the underlying picture portrays a parton 
which is probed by a virtual photon, travels near the light cone, emits a real photon, and remains 
a constituent of the probed hadron. The partonic probability amplitude for such a process is given 
by the generalized parton distribution (GPD) [26, 27, 28]. It has soon been realized that diffractive 
vector meson electroproduction [29], measured in the collider experiments HI and ZEUS at DESY 
(see, e.g., Refs. [30, 31, 32]), can be described in terms of GPDs [33, 34], too. The usefulness of the 
GPDs has also been widely realized in connection with the spin problem, since they encode the 
angular momentum carried by the individual parton species, as explicated by the Ji's sum rule [26]. 
During the last decade they have become an attractive object for theoretical and experimental 
investigations. Today, they are considered as a new concept that provides a link between different 
fields: exclusive and inclusive processes, perturbative and nonperturbative physics (e.g., lattice 
simulations [35, 36, 37, 38, 39, 40] and model building [41, 42, 43, 44, 45, 46, 47, 48, 49]). For 
comprehensive reviews, see Refs. [50, 51]. 

In contrast to parton densities, which depend on the longitudinal momentum fraction of the 
probed parton and on the resolution scale, GPDs encode also transversal degrees of freedom. 
This allows a three-dimensional probabilistic interpretation of the parton distribution either in 
the infinite momentum [52, 53, 54, 55, 56] or in the rest frame [57]. Indeed, to some extent this 
information can already be extracted from present experiments. The main theoretical complication 
arises from the fact that GPDs depend on two longitudinal momentum fractions: x and skewness 
77, which are related to the s- and t-channel exchanges, respectively. The former one is either 
integrated out during the convolution (in the real part of the scattering amplitude) or identified 
with the latter one (in the imaginary part of the scattering amplitude). Even if we precisely knew 
the modulus and the phase of hard exclusive leptoproduction scattering amplitudes 2 , a complete 
reconstruction of the GPD would not be uniquely possible, except with ideal data and employing 
evolution. We note that it was pointed out that already the Fourier transform of the DVCS 

1 This notion, proposed in Ref. [24], stood for expectation values of light-ray operators sandwiched between 
vacuum and final hadronic or between initial and final hadronic states. Today, it denotes the distribution corre- 
sponding to the former expectation value, while the distribution corresponding to the latter one is called generalized 
parton distribution. 

2 We remark that the scan of the GPD shape in a certain region of the momentum fraction plane is possible only 
in the so-called double DVCS process [58, 59, 60] by variation of the incoming photon virtuality and the lepton 
pair mass. 
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amplitude with respect to the skewness parameter provides an image of the target with respect 
to the longitudinal degrees of freedom [61, 62]. The problem of deconvolution can be overcome 
to some extent either by a realistic model, which we at present do not have, or by a hypothesis 
about the form of skewness dependence. In the case of those GPDs that in the forward kinematics 
reduce to parton densities, such a hypothesis can be tested by statistical analysis of inclusive and 
exclusive data, e.g., by means of the x 2 criteria. Certainly, if the skewness problem could be 
solved, the extrapolation to the rj — > case would be simple and would be a step towards the 
experimental access of both Ji's spin sum rule and a three-dimensional picture of the proton. 

Among processes which enable us to access GPDs at present experiments, the DVCS is con- 
sidered the theoretically cleanest one. Indeed, the first experimental DVCS data on the beam 
spin asymmetry in fixed target experiments [63, 64], or the cross section, measured by the HI and 
ZEUS collaborations [65, 66, 67], could be successfully understood even in terms of oversimplified 
GPD ansatze [68]. In contrast, the normalization of the cross sections of vector-meson electro- 
production, predicted to LO in the collinear factorization approach, in general overshoots the HI 
and ZEUS data. This process is widely believed to be affected by power-suppressed contributions. 
They are mainly related to the transversal size of the meson and appear so separately as factorized 
contributions, which might be modelled by themselves [69, 70]. On the other hand, perturbative 
higher-order corrections for this process are large and reduce the size of the predicted cross sections 
[71]. However, to the best of our knowledge, a global NLO analysis of all available experimental 
data has not been achieved so far. In our opinion, this is the only possible way of confronting 
the collinear approach with experimental findings and before this is done, a judgement about this 
approach can hardly be drawn. 

With increasing amount and precision of experimental DVCS data [72, 73, 74, 75], there arises 
a need for a better theoretical understanding of this process. Certainly, for the phenomenology of 
GPDs it is essential to include and estimate perturbative and power-suppressed contributions as 
well as to introduce a flexible parameterization of them. Besides a qualitative or a semianalytic 
understanding of observables in dependence on GPDs, see here, for instance, Refs. [68, 76], also 
fast and stable numerical routines are needed for a fitting procedure. 

So far the perturbative contributions to the DVCS have been worked out to NLO accuracy 
[77, 78, 79, 80, 81], including the evolution [82, 83, 84]. There exist numerical routines in the 
momentum fraction space [85, 86, 87]. In this space the GPD ansatz is given in terms of a spectral 
function [24, 88]. Here one has to model the functional dependencies on two momentum fraction 
variables and the momentum transfer squared. Alternatively, one can work with the conformal 
moments of GPDs [89], which diagonalize the LO evolution equation. This offers the possibility 
of a flexible GPD parameterization, which covers the complete set of degrees of freedom and 
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makes direct contact with lattice results. The problem here is rather to find an appropriate and 
physically motivated truncation in the parameter space. A partial wave expansion with respect to 
the angular momentum, the so-called dual GPD parameterization, has been proposed [90], where 
the partial wave amplitudes are those of the mesonic exchanges in the t-channel with given angular 
momentum. This also allows one to use the angular momentum as an expansion parameter [91] 
and Regge phenomenology as a guideline towards realistic GPD ansatze. 

Unfortunately, if one works with discrete conformal moments, the GPDs are expanded in 
terms of (mathematical) distributions, which live only in the so-called central region (—77 < x < 
77) of the whole support (—1 < x < 1). This is quite analogous to the expansion of parton 
densities with respect to the Dirac function and its derivatives, which live at the point x = 0. 
A variety of approaches have been proposed to resum this formal series and thus to restore the 
correct support of GPDs: smearing the expansion [89, 90], an integral transformation [92, 93], 
Fourier transformation in the light-cone position space [94, 95, 96, 97], and a Sommerfeld- Watson 
transformation [98]. Note that all these approaches should be mathematically equivalent and, in 
particular, those of Refs. [96] and [98] lead to essentially the same representation. So far only 
the smearing method has been extended to NLO accuracy in the MS scheme; however, speed and 
stability remain restricted [99, 100]. 

We believe that the GPD formalism that is based on the Sommerfeld- Wat son transformation 
and Mellin-Barnes integral representation [98] is suitable to satisfy the requirements which we 
spelled out. Moreover, we can employ the power of conformal symmetry to investigate the conver- 
gence of the perturbative series up to NNLO [101, 102]. Therefore, the fivefold goal of this article 
is 

• to present a detailed derivation, based on analyticity and short-distance operator product 
expansion, of the Mellin-Barnes integral representation for the twist-two Compton form 
factors (CFFs), 

• to present the radiative corrections up to NNLO obtained using the predictive power of 
conformal symmetry, as well as to present the radiative corrections in the standard MS 
scheme up to NLO accuracy, including the evolution, 

• to propose an intuitive ansatz for conformal GPD moments, 

• to investigate the convergency properties of the perturbative expansion, 

• to demonstrate the usefulness of the GPD fitting procedure based on this formalism, as well 
as to discuss the requirements for a GPD ansatz. 
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The outline is as follows. In Sect. 2 we introduce the Compton scattering tensor and its decompo- 
sition in leading twist-two Compton form factors (CFFs). Employing their analytic properties, we 
derive dispersion relations for CFFs. In Sect. 3 these dispersion relations together with the short- 
distance operator product expansion are used to represent the CFFs as Mellin-Barnes integrals. 
To diagonalize the evolution operator, we then introduce the conformal partial wave expansion 
of CFFs. In Sect. 4 we employ conformal symmetry in the perturbative QCD sector to find the 
Wilson coefficients for CFFs, expanded up to NNLO in the coupling. We also employ the con- 
formal partial wave expansion for the representation of the NLO corrections in the MS scheme. 
An intuitive ansatz for the conformal GPD moments is introduced in Sect. 5. It is based on the 
internal duality of GPDs and consists of an SO (3) partial wave expansion in the t-channel. In 
Sect. 6 we discuss the convergency of the perturbative series. We then demonstrate in Sect. 7 on 
hand of the DVCS cross section, measured at high energies by the HI and ZEUS collaborations, 
that our formalism is suitable for a fitting procedure that is used to extract GPD information 
from experimental data. Finally, we summarize and conclude. Three appendices contain details 
on our conventions and the evaluation of conformal moments. 



2 Compton scattering tensor and Compton form factors 
2.1 General formalism 

We are interested in the perturbative description of hard photon leptoproduction off a hadronic 
target, e.g., a proton. Besides the Bethe-Heitler bremsstrahlungs process, parameterized by elec- 
tromagnetic form factors, the DVCS process contributes [24, 26, 28]. The amplitude of the latter 
is expressed by the Compton tensor. This tensor is defined in terms of the time-ordered product 
of two electromagnetic currents, sandwiched between the initial and final hadronic state, 

T, u (q,P,A) = ^ Jd A xe^(P 2 ,S 2 \Tj,(x/2)U-x/2)\P 1 ,S 1 ), (1) 

where q = (q± + q 2 )/2 (fi and q 2 refer to the outgoing real photon), while P = P\ + P 2 and 
A = P 2 — Pi. The incoming photon has a large virtuality q\ = — Q 2 and we require that, in the 
limit — q 2 = Q 2 — > oo, the scaling variables 

P ■ q P ■ q 

and the momentum transfer squared A 2 are fixed. The dominant contributions arise then from 
the light-cone singularities of the time-ordered product. This kinematics is a generalization of the 
Bjorken limit, well-known from deep inelastic scattering (DIS). In particular, if the final photon is 
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on-shell, i.e., for the process we are interested in, the skewness parameter r] and the Bjorken-like 
scaling parameter £ are equal to twist-two accuracy, i.e., i] = £ + 0(1/Q 2 ). Note that in the 
following we assume A • P = P| — Pf = 0. 

In the generalized Bjorken limit, we can employ the OPE to evaluate the Compton scattering 
tensor (1). Its dominant part is given by matrix elements of leading twist-two operators. The use 
of the light-cone expansion, a resummed version of the short-distance operator product expansion, 
allows a straightforward evaluation [24] . On the other hand, one might employ the short- distance 
expansion, which yields a Taylor expansion of the Compton scattering tensor with respect to 
the variable l/£. Since this expansion converges only in the unphysical region, i.e., £ > 1, we 
need in addition a dispersion relation that connects the Mellin moments in the physical region 
with the short-distance expansion. This technique well-known from deep inelastic scattering has 
been adopted for nonforward kinematics, e.g., in Ref. [103]. Below we use it within a special 
short-distance OPE, namely, the conformally covariant one. 

Let us introduce a parameterization of the Compton tensor in terms of the so-called Compton 
form factors (CFFs), which has been employed for the evaluation of the differential cross section 
[104, 68]. To leading twist- two accuracy it reads 

T»v(q, P, A) = - ie^qp s 2 + • • • , (3) 



where and e^ vq p = e tlva pq a P^ are the transversal part of the metric 3 and the Levi-Civita 
tensor, respectively, which are contracted X^ u = V^ p X pcr V av with projectors 

a P 

<paP = g af3 _ Q1Q2 / 4 x 
qi ■ <?2 

to ensure current conservation [104]. The ellipsis indicate terms that are finally power suppressed 4 
in the DVCS amplitude, or are determined by the gluon transversity GPD. The latter is a twist-two 
contribution that enters at NLO, evaluated in [105, 81, 106], and is tied to a specific azimuthal 
angular dependence in the cross section [107, 68]. Hence, using an appropriate definition of 
observables, it can be separated from the other twist-two (and twist-three) contributions. For the 
time being we consider its perturbative corrections beyond NLO more as an academic issue and 
will not evaluate it here. 

In the parity even sector the vector 

v° = u(p 2 , s 2 ) (nr + z-^jr) u(Pi, s 1 ) + ---, (5) 



3 The transversal part of the metric tensor reads = — n^n' u — n^n'^, where = (1, 0, 0, ±l)/\/2 and 
n 1 ^ = (1,0, 0, ^fl)/y/2 and at twist- 2 we have n^n^ = q^Puj (q ■ P). 

4 Here we have neglected a further twist-two Compton form factor that can occur for virtual photons that are 
longitudinal polarized. Obviously, it does not contribute to the DVCS amplitude, where the final photon state is 
of course trans versally polarized. 



8 



is decomposed into the target helicity conserving CFF 7i and the helicity flip one 8. Analogously, 
the axial-vector 

A CT = U(P 2 , s 2 ) (n^ l5 + S^f) U(P U 50 + ■ • • , (6) 

is parameterized in terms of 7i and £, where again higher-twist contributions are neglected. The 
normalization of the spinors is U(P, S) , y cr U(P, S) = 2P fT . 
It is convenient to introduce the nomenclature 

with the following choice of arguments: T{y, A 2 , Q 2 ). The variable 

with W 2 = (P 1 + qi )\ u = (P 1 -q 2 ) 2 (8) 



VP 2 £ vP 2 2VP 2 
is in the Breit frame simply given by qtb the conjugate variable of the time Xq. It is an appropriate 
generalization of the photon energy in the rest frame for DIS. Instead of the skewness variable rj, 
depending on u, we will employ the virtual photon mass asymmetry 



= = \ + 0(A 2 /Q 2 ) . (9) 



Obviously, d = in the forward case and d = 1 for DVCS. 



2.2 Analyticity and dispersion relations 

Based on the analyticity of the Compton tensor (1), now we derive a dispersion relation for the 
CFFs, introduced in Eqs. (5) and (6). From the definition of the Compton tensor (1), we read off 
for complex valued q the following properties, see, for instance, Refs. [108, 109], 

[T^(g, P, A; S h S 2 )Y = T Ufl (q*, P, -A; S 2 , S x ) , (10) 
T^q, P, A; Si, S 2 ) = T^(-q, P, A; S u S 2 ) . (11) 

The former of these equations can be viewed as a generalization of the Schwarz reflection principle, 
whereas the latter equation stems from crossing symmetry under the Bose exchange (qi, q 2 , /i, v) <-> 
(— 92, — Qi, //). Furthermore, in the spacelike region, i.e., < Q 2 , and for < —A 2 < — A 2 Iax , 
where the upper limit is constrained by kinematics, the Compton tensor is holomorphic in v except 
for the branch cuts along the real axis. Its absorptive part 

T^(q, P, A; S u S 2 ) - [T u ,(q, P, -A; S 2 , = ^iW^(q, P, A; S u S 2 ) , (12) 
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is given by the commutator of the electromagnetic currents 

W ^=T~A d 4 xe^(P 2 , S 2 \ %{x/2),3 v {-x/2)\ \P X , S,) . (13) 



47re 2 _ 

Finally, using the generalized Schwarz reflection principle (10), we have for the discontinuity, now 
expressed in terms of the more appropriate variable v, 

DiscT^ = T^iv + ie, ■ • • ) - T^iv - ie, • • • ) = AmW^ . (14) 

The spectral representation of the hadronic tensor W^ v and baryon number conservation tell us 



that this discontinuity does not vanish for \v\ > u cut with u cut = (4Q 2 + A 2 )/4VP2, see Fig. 1. 

Obviously, these analytic properties hold for the CFFs, too. The crossing relation (11) implies 
the symmetry relations 

.F v (zv)=^ V HV"-), F A (v,...) = -F A (-v,...). (15) 

To fix the form of the dispersion relation, we should estimate the high-energy behavior of the 
CFFs from the common arguments of Regge phenomenology. The leading meson trajectories give 
rise to a z/ a «( A2 ) behavior with «m(A 2 ) < q;m(0) — 1/2 for A 2 < 0. Thus, for the meson-exchange 
induced part it is appropriate to use Cauchy's theorem with one subtraction, i.e., cf. Fig. 1(a), 



1 r cut ~° dv' 



v 



T(v, 0, A , Q 2 ) = — d> — T{v', 0, A 2 , Q 2 ) + C{&, A 2 , Q 2 ) . (16) 



In the axial-vector case, the subtraction constant C{- ■ ■ ) = T[y = 0, • • ■ ) is zero for symmetry 
reasons, see Eq. (15). In the flavor singlet sector of parity even CFFs T v the pomeron exchange 
appears. It induces a stronger power-like growth with an exponent «p(A 2 ), where «p(A 2 ) < 
ap(0) — 1. For this contribution a second subtraction should in principle be introduced in Eq. 
(16), giving rise to a vC term. However, since a parity even CFFs is a even function in u, this 
constant C is zero. Hence, the analytic properties allow us to derive a singly subtracted dispersion 
relation. Inflating the integration path to infinity, pictured in Fig. 1(b), we are left with a path 
that encircles the cuts on the real axis. Picking up the discontinuity and employing now the 
symmetry (15), we find for all CFFs and |£| > 1 the dispersion relation 5 

, A 2 , Q 2 ) = l - (-^ =f W(£' - tO, 0, A 2 , Q 2 ) + C(0, A 2 , Q 2 ) . (17) 

Here, the upper and lower sign belong to the vector (V) and axial-vector (A) case, respectively, 
where in the latter case the subtraction constant C vanishes. We have also used the fact that the 
imaginary part SmjF(^', • • • ) has the support < 1/(1 + A 2 /4Q 2 ) < 1. 



5 Although the function T is now expressed in terms of £ rather than v, for simplicity we do not introduce a 
new notation. 
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Figure 1: The integration contour used in Eq. (16) and its deformation, employed to derive the 
dispersion relation (17), are sketched in the left and right panel, respectively. Contribution of 
dashed segments vanishes. 

The dispersion relation (17) will be used to relate the short-distance OPE, given as a Taylor 
series with respect to the variable u = l/£ < 1, to the imaginary part of the CFFs in the physical 
region. If we choose to write the Taylor expansion of the CFFs with respect to u as 

oo 

T{t 0, A 2 , Q 2 ) = ^[1 T (-1)'] Ftf, A 2 , Q 2 ) + C(0, A 2 , Q 2 ) , (18) 

3=0 



1 d j+1 1 =F (-IV 

^(■■^STTTlM^aM---) .= n- / ^e^m^-zO,---). (19) 



o 



taking into account the symmetry relation (15), we arrive at 
1 dP +1 

2(j + 1)! du^ 1 ' x ~' "' ' uj=o 2tt 

This relation is valid for < j, where the subtraction term does not contribute. The latter enters 
for j = —1, for which the first term in Eq. (17) drops out, and so we have the trivial identity 

C(tf,A 2 ,Q 2 ) =F v (l/u,$,A 2 ,Q 2 ) . (20) 

u;=0 

3 Operator product expansion to twist-two accuracy 

Using the conformal OPE, presented in Sect. 3.1, and the dispersion relation, we derive in Sect. 3.2 
a Mellin-Barnes representation for CFFs. Then we give a detailed discussion of the subtraction 
constant, which appears in a certain set of CFFs. On hand of a simple toy GPD model we 
investigate some mathematical aspects of the Mellin-Barnes representation. In Sect. 3.3 we present 
the CFFs as conformal partial wave expansion in terms of conformal GPD moments, depending 
on the complex conformal spin. We also provide the rotation of common Mellin GPD moments 
to the conformal ones. 
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3.1 General form of the operator product expansion at twist-two level 



The perturbative QCD predictions for the CFFs (7) might be derived by means of the OPE for the 
time-ordered product of two electromagnetic currents, which is then plugged into the Compton 
tensor (1). The predictive power of conformal symmetry allows for an economical evaluation and is 
the key for the perturbative QCD predictions of the DVCS amplitude beyond NLO, which will be 
presented in Sect. 4. Thus, we choose an operator basis in which this symmetry can be manifestly 
implemented. The evaluation of the CFFs in this so-called conformal OPE is straightforward and 
can be found in Refs. [110, 111, 112, 113, 114, 77, 115]. In this section we mainly present our 
conventions and give a short insight into several aspects of the conformal OPE. 

The time-ordered product of two electromagnetic currents is expanded in the basis of so-called 
conformal operators. The twist-two operators that are built of quark fields read at tree level 



r(3/2)r(i + j; 

2ir(3/2 + j) 



(id. 



y 



1 7+75 J 




(21) 



where, e.g., for three active quarks the nonsinglet and singlet flavor matrices are 



A 



NS 



' 2 


v° 



\ 



and A 



( 1 ^ 
1 
1 j 



f22l 



The gluon operators are expressed by the field strength tensor G a 



fJLV 



duB^ 



gf 



abc r>b 



Gq 

GQ 1 



V 



2J- 2 r(3/2+j) 



id+y- 1 G 




(23) 



The covariant derivatives either in the fundamental or adjoint representation contracted with the 
light-like vector n are denoted as D+=D+ — D+ and d + =d+ + d+ for the total derivative. 
Cj is the Gegenbauer polynomial of order j with index v. The "— " component is obtained by 
contraction with the dual light-like vector n', i.e., n ■ n' — 1. The operators (21) and (23) are the 
ground states of conformal multiplets (towers), which are labelled by the conformal spin j + 2. 
Their Lorentz spin is J = j + 1 and higher states are obtained by multiplying the ground state 
with total derivatives id + . Hence, their spin J is increased by the number of total derivatives. 

The advantage of such a basis is that in the hypothetical conformal limit conformal symmetry 
guarantees that the renormalization procedure can be performed in such a way that only operators 
with the same conformal spin will mix. Beyond LO, this statement is not true in an arbitrary 
scheme, and in particular not in the MS one. Suppose we rely on a scheme in which the conformal 
spin is a good quantum number. Then the only mixing problem occurs in the flavor singlet sector, 
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where quark and gluon operators with the same conformal spin will mix [1 14] . The 2x2 anomalous 
dimension matrix can be simply diagonalized by the transformation 



+/ni \ / W \ 



(24) 



where the 2x2 matrices U l j are expressed in terms of the anomalous dimensions. Here {+, — } 
labels the eigenf unctions of the renormalization group equation in this sector and I G {V, A}. 
Consequently, the renormalization group equation for all operators in question is then simply 

v-jfc&M = -^MaO) ^(aO (25) 

for a G {NS, +, — }. Moreover, if conformal symmetry is implemented in such a manifest way, it 
can be employed to predict the Wilson coefficients of the OPE [110, 111, 112, 113, 114, 77, 115]. 

With respect to the mixing properties of operators, it is convenient to perform also a decompo- 
sition in flavor nonsinglet (NS) and singlet (S) CFFs and express the latter by the two eigenmodes 
a G {+, — } of the evolution operator: 

T = Qls^T + Q 2 S S T , S F = ^ + G f = + (26) 



The fractional charge factors are for three [four] active light quarks 



6 



Q 2 -1 



5 



(27) 



The general form of the short- distance OPE is governed by Lorentz invariance and is given as 
sum over the irreducible representations of local operators. In our conformal operator basis this 
expansion leads to the following leading twist approximation of the CFFs 

F= E ^^Et^^^'H^^j^Q 2 /^^^))^^^,^ 2 ). (28) 

a=NS,± j=0 

Here, they are factorized in the perturbatively calculable Wilson coefficients t , j(t?, Q 2 /fi 2 , cn s (fj,)) 
and nonperturbative reduced matrix elements a Fj(rj, A 2 , fi 2 ), defined below in Eqs. (31) and (32). 
These matrix elements carry conformal spin j + 2 and are nothing but the conformal moments of 
GPDs (see Appendix C). In Eq. (28) the sum runs over the conformal spin and for I=V (A) only 
its odd (even) values contribute (in the following we will lighten our notation by suppressing the 
superscript I G {V, A}). 



6 See Appendix A for more details. Basically, the squared fractional charges Q^s an d Q\ arc obtained from the 
decomposition Q 2 NS NS F + Q| = Q\ U T + Q\ A T + Q 2 S S T [+Q 2 C C T] . 
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The $ = 77/£ dependence of the Wilson coefficients °Cj encodes information on how operators 
with the same spin, but a different number of total derivatives, or different conformal spin projec- 
tion, are arranged 7 . This arrangement is indeed governed in the hypothetical conformal limit by 
conformal symmetry. For instance, to LO approximation this symmetry predicts the form of the 
Wilson coefficients, that are expressed by hypergeometric functions, see e.g., Ref. [114], as 



;i+j)/2,(2 + j)/2 
(5 + 2j)/2 



a cj cc ($) = a cj ee 2 f 1 ( v ^ 1 ; c // i "' v " ; J// " $ 2 ) . (29) 



The normalization a cJ ec is not fixed. However, crucial for our following results is the fact that 
it can be borrowed from deeply inelastic scattering. Namely, setting •d = 0, we realize that a cJ cc 
coincides with the forward Wilson coefficients 



°bf» = a Cf*{$ = 0) = 1 for a = u,d,s,(c). (30) 

Such correspondence is also valid in all orders of perturbation theory, see next section and Ap- 
pendix A. 

The conformal moments Fj = {Hj, Ej, Hj, Ej} arise from the form factor decomposition of the 
matrix elements of conformal operators (21) and (23). The reduced matrix elements 

j-(P 2 ,S 2 \r)]\P 1 ,S 1 ) = U(P 2 ,S 2 )(^H J ( v ,A 2 ,^) 1+ (31) 

j-(P 2 ,S 2 \r)f\P 1 ,S 1 } = U(P 2 ,S 2 )^(r ] ,A 2 ,^ 2 ) 1+ ^ (32) 

of these operators are defined in such a way that in the forward limit we encounter the standard 
Mellin moments of parton densities 

fy(jx 2 ) = lhn a Hj(r), A 2 , ft 2 ) and A a qj {fi 2 ) = lim^^, A 2 , /i 2 ) . (33) 

Note that the matrix elements of conformal operators are tensors of rank j + 1, contracted with 
the light-like vector n. Hence the reduced ones (31) and (32) are polynomials in rj of order j + 1. 
In particular, Hj and Ej are of order j + 1, whereas Hj, Ej and Hj + Ej are of order j, see, e.g., 
Ref. [116]. Moreover, time reversal invariance and hermiticity ensure that all these polynomials 
are even under reflection t] — > —rj. This then implies that for odd j conformal moments, such 
as those appearing in the OPE of JF V form factors, the combination Hj + Ej is actually of order 
J-l- 



7 This can be easily seen from the Taylor expansion of the Wilson coefficients in powers of t]/£. Taking the spin 
J as summation label, the entire expansion is given in terms of inverse powers t;~ J . For each given value J there 
appear J = {1, 2, 3 • • • } matrix elements r] n Fj_ n _i(ri, A 2 , /i 2 ) with n = {0, 1, • • • , J — 1 }. 
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3.2 Employing the dispersion relation 

In this section, we combine the dispersion and the OPE techniques to restore both the imaginary 
and real part of the CFFs in the physical region. We also give an extended discussion of the 
subtraction term. 

3.2.1 Mellin— Barnes representation for CFFs 

We express the Taylor expansion of the CFFs with respect to uj = l/£ in the convenient form 
analogous to that introduced in Eq. (18) 

oo 

A 2 , Q 2 ) = E i 1 ^ ("!)'] W A 2 , Q 2 ) + W A 2 , g 2 ) . (34) 

3=0 

Note that the CFFs and so the Taylor coefficients are physical quantities that are independent of 
the renormalization/factorization scheme that is used for their evaluation. 

First, we evaluate the Taylor coefficients a T^ as defined in Eq. (19), from the conformal OPE 
(COPE) result (28) taking r//£ = rju = d as an independent variable: 

00 1 rl 1 

a ?j * E^i+nto tfV. «.(/*)) $ n ^(A 2 ,/. 2 ) , if = --^F^A 2 ,^) (35) 



n=0 
even 



77=0 



where j is (even) odd in the (axial-)vector case. Here we used the fact that all considered Fj are 

.1 



even polynomials in 77 of the order j or j + 1 and so F- = for I > j or I > j + 1. Consequently, 



the sum on the r.h.s. runs only over even n. 

Next, we use the analytic properties of the CFFs which ensure that the Taylor coefficients (35) 
are given by the Mellin moments of the imaginary part. Corresponding to Eq. (19), they then 
read for general kinematics 

POO 

/ rfa J! 3m^(e,^A 2 ,g 2 ) = 7r^^,A 2 ,g 2 ). (36) 
Jo 

For technical reasons, we extend here and in the following the integration region to infinity using 
Sm"^,---) = for f > 1. 

Now, we can combine Eqs. (35) and (36). For the DVCS kinematics ($ = 1) we have with the 
notation A 2 , Q 2 ) = JF(£, 0=1, A 2 , Q 2 ) 

POO 00 

/ 3mT(^A 2 , Q 2 ) - 7T E ^'+«(27/A ««(^)) a F-+ n (A 2 , /1 2 ) , (37) 

JO „_n 



?i=0 
even 



where we denote the DVCS Wilson coefficients as Cj(Q 2 / fi 2 ,a s (fi)) = Cj(J) = 1, Q 2 /fi 2 , a s (/x)) 
and make use of the photon virtuality Q 2 ~ 2g 2 . It is worth mentioning that the Mellin moments 
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(a) (b) 

Figure 2: (a) The integration contour in Eq. (38) parallel to the imaginary axis in the complex j- 
plane. Adding semicircles results in an integration path which encloses the real axis. Contribution 
of dashed segments vanishes, (b) The integration contour that will be used below in the numerical 
evaluation of CFFs. 

(37) are directly measurable in single spin asymmetries. They contain the sum of all reduced 
matrix elements built by operators with conformal spin equal or larger than j + 2 and so its spin 
is larger than j. 

In the next step we invert the Eq. (37) following the standard mathematical procedure for 
inverse Mellin transform. To do so we analytically continue the l.h.s. with respect to the conformal 
spin, which is trivially done by considering j as complex valued. The Mellin moments are holo- 
morphic functions in j, as long as the ^-integral containing the weight ^(ln^)^ for AT = 0, 1, 2, ... 
exists, and the condition 3?ej > c, with an appropriate constant c, is satisfied. The choice of c 
should exclude the singularities in the complex j-plane arising from the endpoint at £ = 0. Notice 
also that in the limit 3?ej — > oo the Mellin moments of a regular spectral function vanish, except 
for a elastic pole contribution, i.e., a 5(£ — 1) term. The inverse Mellin transform reads 

i rc+ico 

3m<\F(£,tf,A 2 ,Q 2 ) = - / d 3 C J - 1 a ^(tf,A 2 ,Q 2 ). (38) 

Here all singularities of the integrand in the complex j-plane lie to the left of integration contour, 
which is parallel to the imaginary axis, see Fig. 2 (a). The rightmost lying singularity might be 
related to the leading Regge trajectory, and so we have for c = 5te j the condition 8 a(A 2 ) — 1 < c. 

8 The leading Regge trajectory corresponds to £~ Q ( A ") behavior of a T, i.e., £7 -a ( A2 ) behavior of the integrand 
in (37), and the pole in j = a(A 2 ) — 1. Note that in perturbation theory the Wilson coefficients and anomalous 
dimensions possess additional poles at negative integer j. Consequently, if a Tj is perturbatively expanded we have 
the condition max(a(A 2 ) — f , — f ) < c. Moreover, the flavor singlet anomalous dimension matrix in the parity even 
sector has poles at j = 0. These poles could lie right to the pomeron pole j = ap(A 2 ) — 1, which appear in the 
nonperturbative ansatz of GPD moments. Hence, we have in this sector the condition max(ap(A 2 ) — 1, 0) < c. 
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There is one comment in order. The analytic continuation of the integrand in Eq. (38), known 
for non-negative integer j, is not unique. According to Carlson's theorem [117, 118], this ambiguity 
is resolved by specifying the behavior at j ; — >• oo. Consequently, the analytic continuation of the 
Mellin moments must be done in such a way that they vanish (or at least are bounded) in the 
limit Jiej — * oo. Especially, one has to avoid an exponential growth, e.g., by a phase factor 
e ijT_ Obviously, then for £ > 1 the integrand in Eq. (38) does not contribute on the arc with 
infinite radius that surrounds the first and fourth quadrants. Completing the integration path to 
a semicircle, Cauchy's theorem states that the imaginary part of the CFFs vanishes for £ > 1 as 
it should be. 

Finally, to restore the real part of the CFFs from the Mellin moments, we plug the imaginary 
part (38) into the dispersion relation (17), extending the integration interval to < £' < oo, and 
perform the £' integration using the principal value prescription. The existence of the integral is 
ensured by appropriate bounds for dtej. The lower one is max(a(A 2 ) — 1, — 1) < cor max(ap(A 2 ) — 
1,0) < c, see footnote 8, and in addition for the upper one we have c < 1 (c < 0) in the (axial- 
hector case. The integration yields an additional weight factor in the Mellin-Barnes integral (38) 
which is given by tan(7rj/2) and — cot (717/2) for the vector and axial-vector case, respectively. 
The whole amplitude for £ > reads 

^■(tf,A 2 ,Q 2 )+ a C(tf,A 2 ,Q 2 ). (39) 

where the Mellin moments are defined by the series (35). The analogous expression for £ < 
follows from the symmetry relations (15). We stress again that analyticity ties the real and 
imaginary part in a unique way. This is obvious in the axial-vector case, where the subtraction 
constant is zero for symmetry reasons. The vector case will be discussed in the next section. 

By comparing (39) with (34) one sees that with the help of dispersion relation technique the 
analytic continuation of the CFFs in the unphysical region, given by the series (39), yields the 
complex Mellin-Barnes integral (34). The corresponding Mellin moments a Tj are perturbatively 
predicted by another series (35), which depends on the GPD moments a F™ +n . This remaining 
summation we postpone for later. 

3.2.2 Discussing the subtraction constant 

The evaluation of the subtraction constant C from the limit uo — > deserves some additional 
comments. The definition (20) of this constant suggests that it is entirely determined by short- 
distance physics. The relevant local operator has dimension two (or spin zero) and simply does 



<T(^,A 2 ,Q 2 ) = - 



C+lOO 



-J'-l 



i ± 



tan 
cot 



7TJ 

2 
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not exist as a gauge invariant one 9 . So on the first sight one would expect that, as in the forward 
case, the subtraction constant is zero in the (generalized) Bjorken limit. Surprisingly, a more 
careful look at the OPE result (28) shows that all (gauge invariant local) operators contribute, 
because the highest possible order in rj = $/ou, e.g., j + 1, of their conformal moments, cancels the 
suppression factor uji +1 in front of the Wilson coefficients. Hence, the subtraction constant is in 
fact determined by light-cone physics and resumming all terms proportional to Fj 1 ^ leads to: 

oo 

«C(0, A 2 , Q 2 ) ~ 2 Vn-itf, Q 2 /fi 2 , a s ^))r a FS(A 2 , /i 2 ) , °C{$ = 0, A 2 , Q 2 ) ~ . (40) 

n=2 
even 

Since Fj X ' appears only for j + 1 even, this constant vanishes for symmetry reasons in the axial- 
vector case — as it musts — and is absent in the combination 7i + £. Moreover, it is predicted 
to be zero if the virtualities of the photons are equal, i.e., for d = 0. The expression (40) for the 
subtraction constant has been also derived in the momentum fraction representation at tree level 
[119], where it was expressed by a so-called D-term [120]. This term entirely lives in the central 
GPD region and collects the highest possible order terms in t]. Hence, its conformal moments 
provide us if + , e.g., for 77 > 0, 

j\x rfCf (^j D (^j = if +V +1 • (41) 

In the following we discuss the subtraction constant for the CFF £ , rather than 7i. The reason 
for doine; so is that £f +1) originally arises from a nonperturbatively induced helicity-flip, as we 
will see in Sect. 5. We note that such nonperturbative effects are absent in linear combination 
Hj + Ej. 

A subtraction constant is usually viewed as additional information that cannot be obtained 
from the imaginary part. However, if Ej 1 ' is the analytic continuation of Ej, we should indeed 
be able to express the subtraction constant by the imaginary part of the CFF £ . 

Let us demonstrate such a possible cross talk between subtraction constant and imaginary part 
in the most obvious manner. To do so we first make the mathematical assumption that a £(£, $,...) 
vanishes for £ — > 0. In this case the subtraction in the dispersion relation (17) has been overdone. 
Taking the £ — > limit in the dispersion relation indeed leads to the desired relation between 
subtraction constant and imaginary part: 

a C($,...) = - / c^r'Sm^, !?,.■■)■ (42) 

JO 

9 Such an operator might be interpreted as the expectation value of a gauge invariant non-local operator that 
carries spin zero. In the light cone gauge such an operator is for instance given by the local two gluon operator 
g^v B^(Q)B^(0). However, in a gauge invariant scheme its Wilson coefficient is zero. 
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Since the subtraction constant in the equal photon kinematics is vanishing, we derive the sum rule 

a C($ = 0,...) = - / d££~ 1 9tai < E(£,0 = O,...) = O. (43) 
Jo 

We remark, however, that this equality is not of physical interest, because Regge phenomenology 
suggests that the small £ behavior of °€ arises from the pomeron pole and so it will approximately 
grow with l/£ in the high-energy asymptotics. Hence, our mathematical assumption was purely 
academic and the sum rule (43) does not apply. 

Let us now consider a realistic behavior of the CFF 8. Comparing the OPE prediction (40) 
for the subtraction constant with the Mellin moments (35), one realizes that for j — > — 1 the only 
difference is the first term in the series (35). Thus, removing this term and taking the limit leads 
to a formal definition of the subtraction constant 

A 2 ,Q 2 ) ~ 2 lim tejfa A 2 ,Q 2 ) - T7j(tf, Q 2 ///, a s (/i)) a Ef(A 2 ^ 2 )} . (44) 

Here the separate terms on the r.h.s. can be singular and so the analytic continuation plays the 
role of a regularization. We note that, corresponding to the high-energy (£ — > 0) behavior, there 
are also poles and perhaps branch cuts in the complex j -plane, which have to be surrounded first, 
to arrive at Jtej = — 1. This is certainly possible if the Mellin moments are meromorphic functions. 
The situation is more tricky if a Sj contains singularities at j = — 1 and the subtraction procedure 
is maybe not uniquely defined. In particular, the appearance of a fixed singularity at j = —1, e.g., 
a Sj-i proportional contribution, might spoil the relation between the subtraction constant and 
the imaginary part of the corresponding Compton form factor. 

Let us explore the consequences in the case that both scheme dependent quantities t?j=_i and 
a Ej=_ l can be defined by analytic continuation. The GPD moment, appearing on the r.h.s. of Eq. 
(44), might be expressed by the kinematical forward quantity 

a E (0) (A 2 u 2, „ ^ = 0,A 2 ,Q 2 ) 

j { ^ } c^ = o,A2,gv^ a »)' 

see Eq. (35) with & — 0. This provides an important relation, which formally reads 



°C($,A 2 ,Q 2 ) =2 lim 



A 2 Q 2 ) - ^(0, A 2 , QV^ <*■(/*)) a £ . [d = o A 2 Q 2 



(45) 



Remarkably, only the ratio of Wilson coefficients appearing among physical quantities in Eq. (45) 
is potentially plagued by ambiguities. Hence, we must conclude that in the j — > — 1 limit this 
ratio is scale and scheme independent and has the same value in any order of perturbation theory. 
This value can thus be read off from the tree level coefficients (29): 
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This formula has been checked to NLO accuracy. It turns out that radiative a s corrections blow 
up in the limit j — * — 1; however, the leading singularities in j — 1 are d independent. 

Relying on the validity of the prescription (45), we can calculate the subtraction constant from 
the knowledge of the imaginary part. Taking into account Eq. (45), we formally write the desired 
relation as 

°C(0, ...) = - lim ( f°d£ e 3m [<£(£, tf, . . .) - <E(£,0 = 0, . . .)]) , (47) 
1 Uo Jac 

where the integral is analytically regularized. That means that one first has to evaluate the 
integral for max(a(A 2 ) — 1, —1) < Ke j and then perform the analytic continuation to j = —1. 
Certainly, the prescription (47) relies on analyticity. If this property is absent, e.g., because of 
a fixed pole or singularity at j = — 1 with a -^-dependent residue, the subtraction constant a C is 
only partially related to the imaginary part. We remark that in such a case the D-term for GPDs 
[120], introduced as a separate contribution to cure the spectral representation [24, 88], seems to 
be justified on the first view. Otherwise the prescription of Ref. [121] for the restoration of the 
highest possible order in r] is more appropriate 10 . 

Let us remind that the existence of a fixed pole at angular momentum J = 0, i.e., j = 
— 1 was already argued from Regge phenomenology inspired arguments at the end of the sixties 
and it was proposed to access it by the uses of finite energy sum rules [123]. The analyses of 
experimental measurements [124, 125], although not fully conclusive, indicate a fixed pole at 
J = 0. Interestingly, its residue might be related to the Thomson limit of the Compton amplitude 
[126]. On the theoretical side studies in the partonic [127, 128] or light-cone [129, 130] approach 
lead to a real term in the transversal part of the (forward) Compton amplitude that in the language 
of Reggeization was interpreted as a J = fixed pole. Nevertheless, via a subtracted sum rule 
[131, 132, 127, 128] it is related to the imaginary part of the Compton amplitude. This would 
imply that there is no separate .D-term contribution in the central GPD region. 

3.3 Representation of CFFs in terms of conformal moments 

In this section we derive a Mellin-Barnes representation alternative to Eq. (39), which expresses 
the CFFs in terms of conformal GPD moments rather than the usual Mellin moments. The 
advantage of such a representation is that the CFFs are expanded in conformal partial waves, 
which diagonalize the evolution operator to LO and facilitate direct application of conformal 

10 The spectral representation is not uniquely defined [122]. Indeed, if one takes into account that both .D-term 
and spectral function are dependent - they arise from different projections of the same spectral function that enters 
in the improved representation - the choice of a 'common' spectral representation [24, 88] plus D-term is equivalent 
to the improved ones [121]. 
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symmetry beyond LO. In Sect. 3.3.1 we consider the well understood tree-level approximation 
to spell out the mathematical subtleties that appear in the mapping of conformal moments of 
GPDs to the Mellin ^-moments (35) of CFFs. We propose then in Sect. 3.3.2 a Sommerfeld- 
Watson transformation which allows us to derive the desired Mellin-Barnes representation for 
CFFs, which has been already derived previously by other methods [98, 101]. Finally, we discuss 
the Mellin-Barnes representation of conformal GPD moments themselves. 

3.3.1 Lessons from a toy example 

Let us now have a closer look at the series representation (37) of the Mellin moments "jFj, appearing 
in the Mellin-Barnes representation (39) for CFFs. To get some insight into its mathematical prop- 
erties we consider the OPE at tree level. For simplicity, we will discard here the A 2 -dependence, 
which appears only as a dummy variable, as well as the flavor index. The Wilson coefficients of 
the COPE are given in Eq. (29) and for DVCS kinematics, i.e., d = 1, we immediately have 



' V (5 + 2j)/2 j r(3/2)r(3+j; 

If we form partonic matrix elements of the conformal operators, e.g., for the CFF 7i, the conformal 
moments Hj Tee are simply given in terms of Gegenbauer polynomials rfCj /2 (l/r}). According to 
Eq. (35) the corresponding expression of ifj™„ ree reads 

rtree^_r(3/2)r(l + i) , o /2 /l\ , n)tree (-l)"/ 2 r(3/2 + j + n/2)r(l + j + n) 



H tree (n) = K ' ' y ' n 3 C-' - H {n>tree = ^—^ K 1 J 11 v J ' (49) 

1 y n 2'T(3/2 +j) 1 3 \rjj 3+n 2™r(l + j)r(l + n/2)r(3/2 + j + n) ' 1 } 

Consequently, the Mellin moments (35) are for our toy example given by the following series 

w t rec = V (-1) W/2 2 J+1 (3 + 2j + 2n)r(3/2 + j + n/2) 
3 ~ {^(l+J+n)(2+j+n)T(l/2)T(l+j)T(l + n/2Y 1 ) 

even 

As one realizes from the asymptotics of the T-functions, for large n with n> j the summands 
behave as oc n 3 ~ x l 2 and for non-negative integer j this series must be resummed to arrive at a 
meaningful result, which in our case is simply one. Hence, H* ree = 1 and the inverse Mellin trans- 
form yields a (5-function. Note that the exponentially growing factor 2 J in the Wilson coefficients 
(48) is finally cancelled in the conformal Mellin moments (50) by an exponential suppression factor 
in the conformal GPD moments (49). This cancellation ensures that the imaginary part of the 
CFF is vanishing for rj — £ > 1. We conclude that realistic conformal GPD moments for rj > 1 
are exponentially suppressed by a factor 2 _J . 

The real part of the CFF can be obtained from the dispersion relation (17), while the subtrac- 
tion constant follows from (45) 



: tree = lim 2 
i->-i 



ftf* 5 - Cj ree (l)#f )tree =0. (51) 
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This is in agreement with the fact that the highest order terms possible in 77 are absent for our toy 
conformal moments, i.e., H^}^ ee = for n > 0, see Eq. (49). Consequently, the OPE prediction 
(40) shows that the subtraction constant vanishes, too. Thus, we removed the ambiguity caused 
by the subtracted dispersion relation, and arrive at the unique result: 

1 2 
dff J '3mW tree (0 = Ti => 53mft trec (0 = ttS(1 - f) , MeH tiee {£) = - . (52) 



/ 

Jo 



'0 i-e 

Of course, this example is trivial and the findings coincide with the well-known answer, 

W tree (£, # = 1) = f dx ( - 1 , - - ■ 1 . ) H(x, , (53) 



1 



C, — x — ie ^ + x — ie 



arising from the evaluation of the hand-bag diagram with a partonic GPD H(x,rj) = 5(1 — x). 

If evolution is included, an explicit resummation of the series for the Mellin moments (37), 
appearing in the Mellin-Barnes integral (39), cannot be achieved, since each term is labelled by the 
conformal spin and will evolve differently. There are following possibilities to solve this problem: 

• Numerical evaluation of the corresponding series. 

• Relying on an approximation of the series. 

• Arrangement of the integrand in the Mellin-Barnes integral (39) in such a way that the 
integration, instead of Mellin moment, runs over the complex conformal spin 

The third possibility will be worked out in the next section. In the following we shall study 
whether an approximation of the series is appropriate. As we just saw, for complex j with 5te j 
sufficiently small, the £ independent series (49) converges. If the first few terms are numerically 
dominating, we might hope that these terms already induce a good approximation of the CFFs. 
Certainly, if the 77 dependence in the conformal moments is weak in the vicinity of rj — , i.e., if 
the Fj+ n (A 2 ,fi 2 ) are numerically small for larger n, the series might be approximated by a finite 
sum. This should induce a good approximation for the CFFs for smaller values of £. 

Let us shortly demonstrate how this "expansion" works at tree level. We keep in the conformal 
moments (49) only the leading term in if for 77 — ^ 0, i.e., ff^ tTee — \ aric [ a \\ terms with n = 2, 4, • • • 
are neglected. Hence, the Mellin moments (50) read in this approximation 

r«*^®** mm$ - (54) 

The inverse Mellin transform and the use of dispersion relation leads then for £ > to 

>.| 

QfmW tree (0 ~ 9(2 - Q / MeH tree (Z) ~ -2 + C(£ 3/2 ) . (55) 

V 2 - £ 
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The exponential factor 2 J in the approximated moments (54) induces a wrong support. However, 
we have still a useful approximation for small values of £, where the deviation from the correct 
result starts at order £ 3 / 2 . Naively, one might have expected that the accuracy is of order £ 2 . We 
remark, however, that by taking into account the next order in the expansion (49), i.e., n = 2 
as well, the £ 3//2 proportional terms are annulled and the result is then valid up to order (9(£ 5 / 2 ). 
That the expansion of the conformal moments with respect to rj induces a systematic expansion 
of the CFFs in powers of £ is perhaps not true in general. 

3.3.2 Complex conformal partial wave expansion 

We will now change the integration variable of the Mellin-Barnes integral (39) so that the inte- 
gration finally runs over the complex conformal spin and the nonperturbative input is given in 
terms of conformal GPD moments. As mentioned above, in such a basis the evolution operator 
is diagonal with respect to this quantum number. If this diagonality is not preserved in a given 
scheme, we can always perform a scheme transformation so that the conformal GPD moments 
evolve diagonally. Let us start with the Sommerfeld-Watson transformation of the series for the 
Mellin moments Tj into a Mellin-Barnes integral. To do this we rewrite the series 

oo 

V s t& = 1, A 2 ,Q 2 ) = ^(A 2 , Q 2 ) ~ ^2 a C j+n (Q 2 /n 2 ,a s (n)) ^(A 2 ,^ 2 ) (56) 



71 = 

even 



as an contour integral that includes the real positive axis and add two quarters of an infinite circle, 
see Fig. 2, and arrive so at the Mellin-Barnes representation: 

• i-d+ioo i 

^•(A 2 , Q 2 ) ~ - / dn — — a C j+n {Q 2 /fi 2 ,a s (fi)) a F^_ n (A 2 , /i 2 ) , (57) 

where a F™ +n is the analytic continuation of (— l) n / 2a F^ n and the integration contour is chosen so 
that all singularities lie to the left of it. Note that for j = — 1 the expansion coefficient a Fj + ^ is 
absent from this formula. Rather, it is contained in the subtraction constant a C. Plugging this 
representation into Eq. (39), we arrive after the shift of the integration variable j — > j — n at a 
Mellin-Barnes representation for the CFFs 



1 



t(£,a 2 ,q 2 )^- / d 3 r 



c+ioo 



c—ioo 



tan | / Ti j 
- cot | I T 



a C j (Q 2 /fi 2 ,a s (/j,)) a F j (^A 2 ,fi 2 ) . (58) 
Here the conformal moments Fj are now written as a Mellin-Barnes integral 



■ r d+ico ]+ P ij* a F {n) (A 2 u 2 ) 

a F,(£ = r / ,A 2 ,/i 2 ) = - / dn f. , (rj - ie) n J / ]\ 1 , (59) 
JVS 4 J d _ ioo i± e ^-"> w ; sin(ri7r/2) ' v ; 
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valid for < rj < 1. The (lower) upper sign corresponds to the (axial-) vector case. We included 
in this definition a — ie prescription, which ensures that 

iv-ieT = A ^w ^H (60) 



1 ± e*^ - ™) 71 1 ± e i ti~ n ') 7r 

vanishes also for n — > —200 in the complex n-plane rather than generate an exponential j- 
dependent phase factor. We remark that the exponential suppression in the integral (59), arising 
from 1/ sin(n7r/2) in the limit n — > ±ioo, is in general needed to cancel an exponential growing of 

The Mellin-Barnes representation (58) for the CFFs, which we will from now on use in all 
our numerical analysis, has been already derived in two alternative ways in Refs. [98, 101]. As a 
gift of the present approach we found for the conformal moments (31) and (32) the Mellin-Barnes 
integral representation (59), which gives the analytic continuation of the conformal moments. This 
issue has not been completely resolved before. 

We will discuss briefly the representation (59) within our partonic toy GPD to shed light on 
some tricky points. It can be shown that for the conformal moments of GPD H(x,r]) = 5(1 — x) 
the integral (59) exists for < t] < 1, where the integrand can be read off from the second 
expression in Eq. (49). After completing the integration path, it reduces for complex valued j to 
the function ifj ree in Eq. (49) and an addenda that arises from the additional phase factors in the 
definition (59) of conformal moments. For < r/ < 1 we might express our improved findings in 
terms of hypergeometric functions 

= r(3/2)r(3 +i ) + 3 

j v u 2 1+ Jr(3/2 + j) ' V 2 277 J K J 



4 -;-y+2; r(1+j)r(3+j) {3/2 + j/2,2 + j/2 

2^1 



rf - iO 



% + tan(7rj/2) T(3/2 + j)T(5/2 + j) V 5/2 + j 

Hence for positive odd integer j this formula reduces to the Gegenbauer polynomials (times the 
conventional normalization factor). The addenda is proportional to l/(i + tan(7rj/2)), where this 
factor cancels the corresponding one in the conformal partial wave expansion (58). We can again 
close the integration path by an infinite arc so that the first and forth quadrant is encircled. 
It can be shown that the infinite arc does not contribute. Since the integrand stemming from 
the addenda is a holomorphic function in j, we conclude that the addenda finally vanishes for 
< 77 < 1. Consequently, we have found that the appropriate analytic continuation of the 
conformal GPD moments is given by the first term on the r.h.s. of Eq. (61), i.e., the analytic 
continuation of Gegenbauer polynomials is given by the replacement 



r(3/2)r(i+j) m , ^s/2fi\ _ r(3/2)r(3+i) m , ^ f-jj + 3 



2*T(3/2 + j) ^ ° j W ^ 2 1 +^T(3/2 + j) r/J ' Fl V 2 



2rj 



for 77 > . (62) 
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This result allows also for the extension of the r/ region to rj > 1. Of course, symmetry under 
reflection allows for negative values of rj by replacement 77 — > —77. 

3.3.3 Mellin— Barnes representation of conformal GPD moments 

First, notice that a Fj from the Mellin-Barnes representation (59) is not an arbitrary function; 
rather it must be guaranteed that a Fj 

i) is holomorphic in the first and fourth quadrant of the complex j-plane 

ii) reduces for non-negative integers j to a polynomial of order j or j + 1 

iii) possesses for j — > 00 and £ — > 1 an exponential suppression factor 2~ J . 

* (n) 

The requirement i) is satisfied when it holds for a F ■ . To show that also the second requirement 

~(n) 

is fulfilled, let us suppose that the analytic continuation of °F ■ is defined in such a way that it 
does not grow exponentially for n — > 00 with arg(ra) < it/2 and fixed j. Moreover, for simplicity 

a ( n ) 

let us assume that "F^ is a meromorphic function with respect to n. For < rj < 1 we can 
then close the integration contour in Eq. (59) so that the first and forth quadrant are encircled. 
Cauchy's theorem states that a Fj is given as a sum of the residues in these two quadrants 



^AVH^Res 

poles 



1 ± e i(j-n)n s i n ( n7r /2) 



(63) 



This result is in general a series, defining for fixed 77 a holomorphic function with respect to j. For 
non-negative integer j, any contribution that does not arise from the poles at n — 0, 2, 4, • • ■ will 
drop out, due to the factor 1 ± e^ n . Consequently, in this case only the poles on the real n-axis 
contribute and with the requirement a F^ = for n > j + 2 we get the polynomial 

j+p 

a F j (r ] ,A 2 ,fi 2 ) = Y,V na F i ;\A 2 , f i 2 ) for j = {1, 3, 5, • ■ ■ } or j = {0, 2, 4, ■ • • } , (64) 

n=0 
even 

where p = {0, 1}. Note that for p = 1 the highest possible term of order 77 J+1 is now included. 
Moreover, it is required that for rj = £ > 1 the following expression, which appears as integrand 
in the Mellin-Barnes integral (58), 

-i-i 

a F^,A 2 ^ 2 ) ln N j, (65) 



lim 



V<7 V2 



tan 1 / 7rj 
- cot I V T 



vanishes in the limit j — > 00 for arg(j) < tt/2. Here we took into account that the Wilson 
coefficients in fixed order of perturbation theory behave as ln JV j/\fj in the asymptotics we are 
considering. More precisely, l/y/j arises from the factor T(j + 5/2)/T(3 + j) in the normalization 
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of the Wilson coefficients, see Eq. (48), and the logarithmical growing is induced by radiative 
corrections. The condition (65) is necessary to ensure that the support of the imaginary part of 
the CFFs is restricted to |£| < 1. Namely, for £ > 1 we can close the integration path to form a 
contour that encircles the positive real axis. Obviously, the imaginary part drops out and only the 
poles of the trigonometric functions contribute, yielding the series given by the COPE (28). For 
n = £ — > oo we precisely recover the subtraction constant a C, see Eq. (40), which shows that it is 
already contained in the final representation (58) of the CFFs in the conformal moments (59). 

Moreover, we can easily generalize the formula (59) for the computation of the conformal 
moments for a given GPD, which has definite symmetry with respect to the momentum fraction 
x. In this case the analytic continuation of the conformal moments are 



+ f d + io ° l±e ij7T (n-ie) n 2~ n r(3/2 + j - n/2) 



4r(3/2 + j) J d _ loo i ± eti-ny sin ( n7r/2 ) r(i + j - n)r(i + n/2) 

xAC / dxx j - na F(x,n,A 2 ,fi 2 ) 



(66) 



Here the Mellin moments in the second line appear in the form factor decomposition of the 
matrix elements of local operators that are built exclusively out of covariant derivatives. Precisely 
these matrix elements are measured on the lattice [35, 36, 37, 38, 39, 40]. It is understood that 
5te(j — n) is restricted to such values that the integral exists. We note that the polynomiality of 
these moments is ensured by the symmetry we assumed. So for instance for an (anti-) symmetric 
GPD we have for even (odd) values of n 



i nl n\ n +P 

- dxx na F{x,n,A 2 ^ 2 ) = dxx na F(x,n,A 2 ^ 2 ) = y2 a F m (A 2 ,^ 

1 J-l Jo i=0 



n\ (67) 



where p takes the values and 1 for even and odd n, respectively. However, irrespective of the 
symmetry p = holds true for the combination H + E, H, and E. 



4 DVCS amplitude up to next-to-next-to-leading order 

In Sect. 4.1 we first discuss the predictive power of conformal symmetry in a conformally invariant 
theory for general kinematics. Here we rely on a special subtraction scheme, the so-called conformal 
one, in which the renormalized operators are covariant under collinear conformal transformations. 
This allows us to predict the functional form of the Wilson coefficients up to a normalization 
factor. This factor is simply the Wilson coefficient that appears in DIS. To reduce the massless 
QCD to a conformally invariant theory, we intermediately assume that the QCD has a non-trivial 
fixed-point so that the trace anomaly vanishes. Then we move to the real world and discuss the 
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inclusion of the running coupling. In Sect. 4.2 we present the results for the DVCS kinematics. 
The recent progress in the evaluation of radiative corrections in DIS allows us to present the CFFs 
to NNLO accuracy. However, the trace anomaly induces a mixing of conformal GPD moments due 
to the evolution. Also note that the forward anomalous dimensions of parity odd conformal GPD 
moments are mostly unknown beyond NLO. In Sect. 4.3 we present the NLO corrections to the 
twist-two CFFs in the MS scheme and discuss the inclusion of mixing effects in the Mellin-Barnes 
representation. 



4.1 Conformal OPE prediction and beyond 

In the hypothetical conformal limit, i.e., if the (3 function of QCD has a non-trivial fixed point 
a* 7^ 0, the 7//£ dependence of the Wilson coefficients is predicted by conformal symmetry [110, 
111, 112, 113]. For all the cases we are considering, Wilson coefficients a C I , with a = {NS, +, — } 
and I = {V, A}, have the same general form. Thus, when there is no possibility of confusion, we 
will simplify our notation by suppressing both of these superscripts. The COPE prediction for 
general kinematics reads [114] 



-(2 + 2j + 7j (a;))/4,(4+2j + 7;f (a;))/4 r,-\ fn 



w(,w,«o^(a.vf ■ - ■ ;5?2;;;> : ))/2 uv ~*""$ SP*" 2 - (68 » 

where the normalization Cj remains so far unknown. In the kinematical forward limit (r/ = 0) the 
Wilson coefficients coincide with those known from DIS 

7j«)/2 



(69) 



see Appendix A. Neither Wilson coefficients (68) nor conformal GPD moments (31), (32) mix 
under collinear conformal transformations. This implies that they evolve autonomously with 
respect to a scale change, 

^j-F 3 {...^ 2 ) = -tKIF^..,^ 2 ), (70) 

/i|-Q(...,g 2 //i 2 ,a:)] = C i („.,Q a /Ai 2 J a:)7i(a:). (71) 

Here the anomalous dimensions are the same as in DIS. We also realize that the conformal spin 
j + 2 is a good quantum number. 

Unfortunately, conformal symmetry is not manifest in a general factorization scheme, and in 
particular not in the MS scheme. What is required is a special scheme, which ensures that the 
renormalized conformal operators form an irreducible representation of the collinear conformal 
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group SO (2,1) [114]. The breaking of the covariant behavior with respect to dilatation, e.g., shows 
up in the mixing of operators' expectation values and of Wilson coefficients under renormalization: 

7 j 

A4 _F i (...,/* 2 ) = -J2ljd<Xs(ri)v j - k Fk(-,V 2 ), (72) 



da 

r k=0 



V-C 3 (...,Q 2 /v 2 ,a s (f,)) = ^C h (...,(?/f,aM)^(aM) (jj , (73) 

k=j 

respectively 11 . This mixing arises in NLO and is induced by the breaking of special conformal 
symmetry in LO. The mixing matrix, i.e., 7^ for j > k, contains besides a (3 proportional term 
also one that does not vanish in the hypothetical conformal limit. Knowing the form of the 
special conformal symmetry breaking in the MS scheme to LO, it was shown that, after rotation, 
the diagrammatic evaluation of NLO corrections to (flavor nonsinglet) anomalous dimensions 
[133, 134, 135] and Wilson coefficients [78, 79, 80] in the MS scheme coincide, for (3 = 0, with the 
conformal symmetry predictions [136, 114, 77]. A formal proof about the restoration of conformal 
symmetry, valid to all orders of perturbation theory, has been given in Ref. [114]. Hence, we 
have no doubt that conformal symmetry can be effectively used for the evaluation of higher order 
perturbative corrections [137, 101, 102]. The interested reader can find a comprehensive review of 
the uses of conformal symmetry in QCD in Ref. [115]. 

Let us suppose that we employ a scheme in which conformal symmetry is manifest in the 
hypothetical conformal limit. Such a scheme we call a conformal subtraction (CS) scheme. We 
discuss now the structure of the conformal OPE beyond this limit, i.e., the implementation of 
the running coupling constant in the Wilson coefficients (68) and the evolution equation. In 
the simplest case, i.e., forward kinematics, the conformal representation is trivial, as there are 
only operators without total derivatives left, so we can restore the Wilson coefficients of the full 
perturbative theory from the result (69) of the conformal limit. Namely, renormalization group 
invariance of the CFFs allows us to revive the /3-proportional term by expanding the exponential 12 

CiK) (^) 7jKV2 =► cjMQ)) exp j jT ^7iM/0)} (74) 

in terms of ln(Q 2 /fi 2 ). Of course, also the evolution of Mellin moments of parton densities is 
governed by a diagonal evolution equation, i.e., by Eq. (72) with rj = 0. Finally, we are using the 



11 Both these equations together ensure that the CFF^£.£ J 1 Cj(- ■ ■ ,Q 2 /fj, 2 ,a s (iJ,))Fj(- ■ ■ , fj,) is a renormaliza- 
tion group invariant quantity. 

12 We consider Cj and jj as (3 independent quantities and so their perturbative expansion in the fixed-point 
and full theory looks the same. Certainly, both of them contain to two-loop accuracy (3q proportional terms and 
the anomalous dimensions to three-loop accuracy might be also rewritten in terms of the (3 function expanded to 
two-loop accuracy. 
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normalization condition 

c J (a s )=cf S (a s ) (75) 

and thus we end up with the standard result for the Wilson coefficients of DIS, evaluated in the 
MS scheme. 

For general kinematics the /3 proportional term is not automatically fixed, so different pre- 
scriptions might be used. For a discussion of this issue see Ref. [137]. We will point out here two 
appealing possibilities. 

Let us first assume that the conformal symmetry breaking can be entirely incorporated in the 
running of the coupling so that the evolution equation for conformal GPD moments is diagonal: 

^Fjirj, A 2 , /i 2 ) = - 7 > s (/i)) Fj ( V , A 2 , /x 2 ) . (76) 

Then also the Wilson coefficients in the full theory are simply obtained by replacing the expression 
for the resummed evolution logs. However, it is likely that the trace anomaly generates a non 
vanishing (3 proportional addenda, which will appear at two-loop level. If it does, renormalization 
group invariance fixes the scale dependence and so the Wilson coefficients in such a scheme read 
for the full theory 



CMS, 1, a s (Q)) + ^( as (g))AC i (7 7 /e, 1, 

where Cj(r)/£, l,a s ) is the COPE prediction (68). In the forward limit we again require that the 
Wilson coefficients coincide with the DIS ones, evaluated in the MS. Thus, we have the condition 
that the (3 proportional addenda vanishes, i.e., 

limAC,(r / /e,l,a s (Q)) = 0. (78) 

The disadvantage of such a conformal scheme is that a consistent determination of the f3 pro- 
portional addenda ACj requires a diagrammatical evaluation at NNLO, e.g., in the MS scheme. 
Both the knowledge of the /^-proportional terms in the two-loop corrections to the Wilson co- 
efficients and the off-diagonal three-loop corrections to the anomalous dimensions are required. 
While the former are known in the MS scheme [138, 139], the calculation of the latter, although 
reducible to a two-loop problem by means of conformal constraints, still requires a great effort. 
Fortunately, to NNLO only the quark sector suffers from this uncertainty. The best what we 
can do to NNLO accuracy is to minimize the effects of this lack of knowledge. To do so, we will 
perform a transformation into a new scheme, called CS, where the normalization condition reads 

AC j { V /t,l,a a {Q)) = 0, (79) 
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while the mixing term is now shifted to the evolution equation. This mixing term will be suppressed 
at the input scale by an appropriate initial condition and we might expect that the evolution effects 
are small, see discussion in Sect. 6.3 below. Hence the Wilson coefficients read now 

oo 

Cjiv/t <W. <*.Qi)) = C k (rj/^ 1 , a s (Q)) (80) 

k=j 



xPexp 




nfj(ai a (n'))6 k j + ( -J -A kj (a s (fi')) 



H^M ■ ■ ,f) = -7i(a.(A*)) F j (■ ■ ■ ,/x 2 ) - ^^^^-*A iJb (a.(Ai))F fc (. • • , fi 2 ) . (81) 



where Ck(r)/£,l,a s (Q)) is given in Eq. (68), A k j is determined by the off-diagonal part of the 
anomalous dimension matrix, and V denotes the path ordering operation 13 . The renormalization 
group equation for the conformal moments reads 

Since this mixing term appears for the first time at NNLO accuracy, and is suppressed at the 
input scale by the initial condition, we might expect that its effects will be small. 

Formulae (80) and (81) are our main result that can be used for the numerical investigation 
of radiative corrections to DVCS up to NNLO. 

4.2 Perturbative expansion of Compton form factors 

For DVCS kinematics the general expression for the Wilson coefficients (80) simplifies considerably. 
Since to twist-two accuracy 'd = t]/^ = 1 is valid, the argument of hypergeometric functions is 
unity and they simplify to products of V functions. Having also in mind that in DVCS the photon 
virtuality Q 2 ~ 2Q 2 is considered as the relevant scale, the COPE prediction (68) together with 
Eq. (80) leads to the following expression for the DVCS Wilson coefficients 

Q(Q 2 //i 2 , a a (ji))) = C0 = 1, Q 2 /(2/i 2 ), a B (fi)) (82) 

= Cj(a s (Q)) W 3\ F /o, •_, i / rrm \ exp i / "7X7? ( a s(/^ )) \ 

r(f) T(3 + j + 27j(<MQ))) Us V ) 

at the considered NNLO accuracy. Note that the off-diagonal part of the evolution operator does 
not appear in this approximation, and that we have slightly changed the normalization condition 



13 



The path ordering operation defined by 



Pexpj r ^L ljk (a s (iJ.'))\ = 1+ ^r7 jk (a s (fJ,')) + r% ^ %5Z7im(a.0*'))7mfc(a.(A*")) + --- 

kJQ^ 1 ) J Q ^ J Q /* J Q rn—k 

enables the compact representation of the solution of the non-diagonal renormalization group equations. Obviously, 
in diagonal case, i.e. for 7^ = Sjklj, path ordering converts the evolution operator in a simple exponential function. 
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for fi = Q. In fact to arrive at the above result, Eq. (80) has to be multiplied by a factor 

(V2)*- (a » exp {- jf^ ^(,.00)} = 1 + If §7f ln'(2) + 0(a]) , (83) 

and such terms can be compensated by a change of the normalization condition (79). Here and in 
the following the first expansion coefficient of (3(g) /g = (a s /4:7r)(3 + 0(a 2 s ) is (3 = (2/3) — 11, 
where n/ is the number of active quarks. 

We now use perturbation theory to determine Cj(Q 2 / /i 2 , cn s (/f)) to NNLO accuracy. In contrast 
to Fj(- ■ ■ ,/i 2 ), it is customary to express Cj(Q 2 / (i 2 ,a s (fi)) completely expanded in a s (fi), i.e., 
without resummation of leading logs from /x to Q. This actually corresponds to the result that 
one would obtain from a diagrammatical calculation, without using the renormalization group 
equation. Obviously, since the leading logarithms in Fj will be resumed and in Cj not, one will 
end up with a residual dependence on the factorization scale \x. This scale might also be used 
as the relevant scale in the running coupling constant. However, for generality, another scale /x r 
is better suited for the expansion parameter a s (fi r ). The truncation of perturbation theory then 
implies that our results will depend also on this renormalization scale /i r . 

The perturbative expansion of the DVCS Wilson coefficients (82) in terms of a s (/j r ) is done in 
a straightforward manner by means of 



a 2 (u ) u 2 

a s (ii) = a s (fx r ) + Po^f 11 In ^ (84) 



and consequently we can expand the exponential factor 

2 \ 7j(as(Mr))/2 |- 2/., \o r> 2 /, 2 D 2 

(2tt) z 8 J fi z \i% 



Furthermore, we use the expansion 



a^rd + j + H-) (» 2 Y /2 _ 2^r(5/2 + j) ^ s^\q 2 /^ 2 ) 



r(|)r(3 + j + | 7j ) \Q 2 J r(3/2)r(3 + i) ^ 2-m! 

where the so-called shift coefficients, which also include the factorization logs, are defined as 

^(o 2 /n 2 ) - £L ^ 2 V ^r(3 + j)r(5/2 + J + p) 

j U/M dff»\&) r(5/2 + j)r(3 + j + p) p=o' [ * n 
The first two shift coefficients read in terms of harmonic sums as 

sf ) (Q 2 /fi 2 ) = S 1 (j + 3/2)-S 1 (j + 2) + 2\n(2)-\n%, (88a) 

sf(Q 2 /fi 2 ) = (sf(Q 2 /f)f - S 2 (j + 3/2) + S 2 (j + 2), (88b) 
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where the analytical continuation of these sums are defined by 



S 1 (z) 
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\nT(z + 1) + 7^ and S 2 (z) 
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9) 



dz dz 2 
with -f E = 0.5772 . . . being the Euler constant. Taking into account the perturbative expansion 
of anomalous dimensions and DIS Wilson coefficients, written as 
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the DVCS Wilson coefficients (82) take the form 
T+ l V{j + 5/2) 



where 
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Notice that the renormalization group logs ln(Q 2 //i 2 ), caused by the running of a s , and propor- 
tional to po, appear for the first time in the NNLO correction C- . 

4.2.1 Flavor nonsinglet Wilson coefficients 

The flavor nonsinglet Wilson coefficients have the form as indicated in Eqs. (91a)-(91c). We only 
have to decorate them and the anomalous dimensions with the corresponding superscripts. It has 
been already taken into account that the DIS Wilson coefficients are normalized as 



NS 1(0) 



(92) 



The radiative corrections to NLO involve the Wilson coefficients of the DIS unpolarized structure 
function F%, 



NS V(l) 



Sii 1 +3) + n S iU + 2 ) - o + 



2 "" ' 2 ' 2(j + l)(j + 2) 
for the vector case and those of the polarized structure function gi, 



NS A(l) 
C 3 



S 2 (l+j) + h 1 (j + 2)-^- + 



3-2ff 1 (j) 
2(j + 1)(j + 2) 



S2U + K 



(93) 



(94) 
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for the axial-vector one. The LO anomalous dimensions read in both cases 



NS 1(0) 



'77' =-Cf[S+ g - 1} 2 (j - 2) - 45x0- + 1) ) , (95) 



where Cjr = (Af? — l)/(2A r c ) with iV c being the number of colors. 

The expressions for two-loop quantities are lengthy. The nonsinglet anomalous dimensions 

for the vector and axial- vector case are given, e.g., in Ref. [140], by analytic continuation of j-odd 

(2) 

and j-even expressions, respectively. The NNLO Wilson coefficients Cj can be read off for the 
vector and axial- vector cases from Refs. [141] and [142], respectively. Implementation of evolution 
to NNLO accuracy requires also the anomalous dimensions to three-loop order [143]. Instead of 
the exact expressions one might alternatively use an analytic approximation, given for the vector 
quantities in Refs. [144, 143]. 

4.2.2 Flavor singlet Wilson coefficients 

So far the {+, — } basis, which resolves the mixing problem in the flavor singlet sector, has been 
used for presenting the conformal predictions. However, the Wilson coefficients are in the literature 
usually given in the {Q = E, G} basis. For convenience, we might express the flavor singlet results 
in the basis of quark and gluon degrees of freedom. We write the Wilson coefficients as a two- 
dimensional row vector 

C;.(Q 2 //i 2 , «.(a0) = , G Cj) O.(a0) for I = {V, A} . (96) 

In the following we drop the superscript I. The rotation between the two bases is given by the 
transformation matrix Uj, defined in Eq. (24). Hence, for the Wilson coefficients we have 

C^QV/i 2 , o^)) = ( + Q , -C s ) (Q 2 //i 2 , a.(Ai)) Ufa, Q) , (97) 

where the {+, — } entries are given in Eq. (82). Using the property (B.4) of the transformation 
matrix within /xq = Q, derived in Appendix B, allows us to perform the transformation at the 
normalization point \x = Q and afterwards restore the scale dependence by an backward evolution: 

C j (Q 2 /»\a 8 (ji))=( + C j , -C J )(Q 2 /fi 2 = l,a s (Q))U J (Q,Q)£j l ( f i,Q). (98) 

Here E~ x is the inverse of the evolution operator in the {E, G} basis, defined below in Eq. (120). 
Since Uj(Q, Q) rotates the anomalous dimensions [depending on a s (Q)], we can express for \i—Q, 
the flavor singlet contributions in the form 

c,«JV = WO) = rp/Sstl) l ^"'^'" '' H* ■ «* = • °^ • 
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where the anomalous dimension matrix is defined as 



7, 



2^ 



E5^(0) SG 7 (0) 

GeJo) G« (0) 
' 3 1 3 



0(al 



(100) 



This matrix valued expansion, although resulting from somewhat more involved derivation, is in 
complete analogy to the flavor nonsinglet result (86). 

Acting with the evolution operator £j 1 (/x, Q) on Eq. (99) removes the constraint Q 2 /fi 2 = 1 
of the shift functions argument. As in the preceding section, we consequently expand with respect 
to a s (fi r ) and finally write the Wilson coefficients (96) in complete analogy to Eq. (91) as 



T +l T(j + 5/2) 
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To LO the DIS Wilson coefficients are normalized as 
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In NLO the quark Wilson coefficients correspond to the non-singlet ones given in Eqs. (93) and 
(94), 
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while the gluon ones read 
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The entries of the anomalous dimension matrices read to LO in the vector case: 
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and in the axial- vector case: 
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where Ca = N c and Tp — 1/2. At this order the anomalous dimensions in the quark-quark 
channels are identical to the flavor nonsinglet ones 

(113) 

given in Eq. (95). The NNLO DIS Wilson coefficients for the structure functions F\ and g\ are 
given in Refs. [141, 142]. The NLO singlet anomalous dimensions for the vector and axial- vector 
cases can be found in Refs. [140, 145] and [146, 147], respectively. To treat the evolution to NNLO 
approximation we need also the NNLO flavor singlet anomalous dimension matrix. It has been 
calculated for the vector case in Ref. [148]. We remark that all these quantities are expressed in 
terms of rational functions and harmonic sums. Several numerical routines for them are available, 
e.g., in Refs. [149, 150, 151, 152]. Here again one might rely on analytic approximations for the 
quantities in question, see Refs. [153, 148]. 



4.2.3 Expansion of the evolution operator 

The evolution of the flavor nonsinglet (integer) conformal moments in this CS scheme is governed 
by 
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The solution of the renormalization group equation is given by the path-ordered exponential, 
appearing in Eq. (80). Unfortunately, the mixing matrix A^ fc s is not completely known, and so we 
can here only deal with the solution for A^. = 0: 

F j(£, A 2 , /i 2 ) = £,•(//, /io)-Fj(£, A 2 , nl) , where £,•(//, fi ) = exp ^ - / —jj{fi') > . (115) 

In the numerical analysis we will resum only the leading logarithms and expand the non-leading 
ones. The result, see Ref. [137], reads 
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The expansion coefficients of the (3 function are 
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The flavor singlet case can be treated in the diagonal {+, — } basis in the analogous way as 
the nonsinglet one. Here we present the solution in the {Q, G} basis: 
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We will expand the evolution operator 
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up to order a 2 , and while leading logs remain resummed the non- leading ones are expanded. To 
have a condensed notation, we perturbatively expand the result in terms of the leading order 
{+, — } modes: 
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Here the projectors on the {+, — } modes are 
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where the eigenvalues of the LO anomalous dimension matrix (i.e., from Eq. (B.2)) are 
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*Here the path ordering V is related to the non-diagonality of the anomalous matrix in {Q, G} basis. 
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A straightforward calculation leads to the matrix valued coefficients 
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where the /i dependence is accumulated in the functions ab Rj(/i, Ho\ri) = ab Rjj(fi, [i \n) defined by 
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The expansion of the evolution operator (116) or (121) will then be consistently combined with 
the Wilson coefficients (91) or (101), respectively, see for instance Ref. [137]. 



4.3 MS results up to NLO order 

The DVCS radiative corrections in the MS scheme are known only to NLO. The NLO corrections 
to the hard-scattering amplitude have been obtained by rotation from the conformal prediction 
[114, 77] and agree with the diagrammatical evaluation of Refs. [78, 79, 80]. Employing Eq. (C.6), 
the corresponding Wilson coefficients can be obtained by determining the conformal moments of 
the hard-scattering amplitude, e.g., summarized in Ref. [154]. 

The integrals which are needed for the quark part are given in Appendix C of Ref. [137] for 
integer conformal spin (see also the last paragraph in Appendix C.l). The analytical continuation 
is straightforward and so in the MS scheme we have 



E C/ (1) (Q/^ 



2S 2 (l+j) 



2S?(l+£ 



9 , 5-45' 1 (j + l) 



- + 



1 



2 2(j + l)(j + 2) ( J + 1)2( J + 2 ) 2 J 
9 3-45i0' + l 



+ 



2 Q 2 



(127) 



+ 



1 



2 2(j + 1)(j + 2) ( J + l)2( J + 2 ) 



S V 0) a 2 



'1281 



in vector and axial- vector case, respectively, while N ^C^ 1 ' ) = s Cj ■ 

The integral expressions needed for the determination of the singlet gluon contributions we 
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list in App. C.2. For the gluon conformal moments we obtain 
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The complete anomalous dimension matrix in MS scheme is known to two-loop accuracy 
[136, 155, 82]. We present here the more involved case of the evolution in singlet sector and will 
just state the substitutions needed to obtain nonsinglet case results; for the extensive account of 
the singlet case evolution see also Ref. [100]. Note that our normalization, adopted from DIS, 
differs from the original one. Anomalous dimensions used in this work are related to those in Ref. 
[100], denoted as 

~r jk = N ]7jk \ [m N- k 1 . (131) 

The transformation matrix is 

Nj = N U) ( o 6 ) ' (132) 

with the normalization factor 
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originating from the altered definitions (21) and (23) of conformal operators. Explicitly, we have 
for non-diagonal entries 
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coincide with the DIS anomalous dimensions. 

The evolution operator in the MS-scheme leads already at NLO to a mixing of conformal GPD 
moments. For integer conformal spin these moments evolve as 

Fj(ri, A 2 , fi) = ^2 1 T ( 2 _1) S jk (fi, /i ; v)F k (v, A 2 , //<,) • (136) 

k=0 
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To NLO accuracy the evolution operator is expanded as 15 
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The diagonal term £jj(li, i^o] rj) coincides with the evolution operator £j(ji,ji ) of the preceding 
section, where the diagonal anomalous dimensions 7 • are the same in MS and CS scheme. Also 
the projector a Pj is the same as before and is defined in Eq. (122). The matrix a $B^(/i, /io), 
defined for j — k = {2, 4, ■ • • } (while otherwise), assumes the form analogous to (124) 

ah B§ = - a % k (fi, /i | 1) a P 3 7$ b Pk , for j - k even , (138) 

while 7^ can be expressed as 

^=h {0 \9 + d((3o-7 m )] jk . (139) 

This equation holds for parity even and odd anomalous dimension matrices, where g^ k and dj k 
transform according to Eq. (131). In analogy to Eq. (134), the matrices gj k , and dj k can be read 
off from the results given in Ref. [100]. Taking into account the usual properties of projectors, one 
can conveniently rewrite ^^(fi, fi ) as 

a ^Wo) = - a ^M/^o|l) {%-%) [(A) - %) a Pj d jk b P k + a P 3 g jk b P k ] . (140) 

The NLO solution of the evolution operator convoluted with the Wilson coefficients can be more 
easily treated as the evolution of the conformal GPD moments itself. Namely, the convolution, 
represented as an infinite series and understood as an analytic function of j, can be expressed in 
two equivalent forms 
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15 If one would be interested to resum all logs rather than only the leading ones, one certainly has to deal also 
with the rcsummation of the non-diagonal entries. This problem has been solved in the flavor nonsinglct sector for 



r) = 1 in Ref. [156]. 
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The form on the r.h.s. corresponds to "shifting" the evolution to Wilson coefficient which is more 
convenient for numerical treatment of non-diagonal terms. In this representation it is also obvious 
that in DVCS kinematics the mixing under evolution is not suppressed by powers of £ 2 . The 
evolved Wilson coefficients can furthermore be decomposed as 

oo 

even 

(142) 

Of course, we understand that a consequent expansion in a s to NLO will be performed. To have 
a numerically more efficient treatment this sum can be transformed into a Mellin-Barnes integral. 
Employing the fact that the only residue of cot(7r/c/2) is 2/tt for even integer values of k, we 
straightforwardly arrive at 

C 3 (Q 2 /fi 2 ^^ ) = C j (Q 2 /fi 2 ,a.(jj,))e jj (ji,fi ;l) (143) 

1 f c+ico ( nk\ 

-— J dk cot ( — J C j+k+2 (Q 2 /fi 2 ,a s (fi)) £ j+k+2 ,j(fi, fi ; 1) , 

where —2 < c < 0. 

The formulae for the flavor nonsinglet sector (see, for example, Ref. [137]) are simply obtained 
by reducing the matrix valued quantities to real (or complex) valued ones, by performing the 
replacements 

■y . — > N ^-y . «\ ■ -> NS -Y (0) = ES -Y (0) a P ■ -> 1 

' 3 '3 ' A 3 1 3 1 3 ' 3 ^ 1 ' 

7j* Ijk = Ijk 7 9jk ~* 9jk = 9jk , djk -> rfjfc = rfj-fc . (144) 

We recall that in our convention the off-diagonal entries (j > k) are dressed with a factor Nj/N k , 
e.g., NS 7jfc = (Nj / j k \^iQQ^, and that only to NLO accuracy, considered here, NS 7jfc coincides 
with the quark-quark entry EI 7jfc- 



5 Parameterization of GPDs 

For the phenomenology of GPDs it is crucial that one has a realistic ansatz at hand. At present, it 
is popular to use a GPD ansatz [33] that is based on the spectral representation of GPDs [24, 28] , 
for a review on this subject see [157]. We think that our formalism offers the possibility to approach 
this problem from a new perspective and that it finally leads to a more flexible parameterization. 
It is also remarkable that the task of modelling the x dependence of GPDs by using the constraints 
of the lowest moment [158, 159, 160], i.e., j = yields elastic form factors, turns in our approach 



40 



into the inverse problem. Namely, it is rather the (conformal) spin dependence of the form factors 
that determines the x dependence, where j = serves as an 'boundary' condition. Since it is 
timely to find a parameterization that is the most appropriate for a fitting procedure, we here 
only briefly outline our considerations. 

We recall that the GPD support can be separated into two different momentum fraction regions, 
both having a dual partonic interpretation, namely, as the exchange of partons in the s-channel 
and as mesonlike configurations in the t-channel. Moreover, these regions are tied owing to the 
Poincare invariance, see, e.g., Ref. [98]. This well-understood support property and its partonic 
interpretation might be considered a reflection of the conjectured duality between t- and s-channel 
processes; for a review of its phenomenological application, see Ref. [161]. This point of view opens 
the door to employ strong interaction phenomenology (combined with some basic principles) in 
finding a realistic model ansatz for GPDs and partial extraction of its parameters. We believe 
that it is a promising starting point to think in terms of t-channel exchanges combined with SO (3) 
partial wave analysis and Regge phenomenology and then to perform the crossing to the DVCS 
process. This approach is complementary to those which are based on modelling GPDs in terms 
of partonic degrees of freedom, see, for instance, [41, 42, 43, 44, 45, 46, 47, 48, 49]. In the next two 
sections we outline the SO (3) partial wave analysis of DVCS and, based on an intuitive picture, 
propose an ansatz for the partial wave amplitudes, appearing in conformal GPD moments. In the 
third section we try to get insight into the properties of this ansatz, and, finally, in Eq. (165) we 
present its approximated version convenient for numerical studies and used in subsequent sections. 
The reader more interested in the applications might just skip to that point. 

5.1 SO (3) partial wave decomposition of conformal GPD moments 

It has already been proposed in Ref. [90] to decompose conformal moments of meson GPDs into 
irreducible SO(3) representations, namely, in terms of Legendre polynomials labelled by the orbital 
angular momentum quantum number. These functions enter the partial wave decomposition of 
the t-channel scattering amplitude and depend on the scattering angle 9 = 6 cm , defined in the 
center-of-mass frame, see Fig. 3(a). After crossing, cos# becomes 16 the inverse of the scaling 
variable —77, i.e., cos 6 = — 1/rj [24, 90]. The basis of Legendre polynomials has later been adopted 
for the conformal moments of nucleon GPDs, too [91]. Let us here start with the analysis and 
discussion which basis of polynomials one should employ. 

To avoid cumbersome technical details, we restrict ourselves here to massless hadrons, so that 
the crossing relations between helicity amplitudes become simple: the signs of both four-momenta 

16 For the sake of simplicity, here we neglect \J 1 + M 2 /P 2 proportional corrections and those which die out in 
the Bjorken limit. 
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Figure 3: The fusion of two photons into a pair of (massless) mesons in the center-of-mass frame 
(a), partonic space-time picture (b), and hadronic point of view (c). The thick arrows indicate 
the spin projection of the particles. 

and helicities of the crossed particles are to be changed. However, one has to be careful with 
the virtual photon 17 , which we will not cross, but whose virtuality causes a rotation of helicity 
amplitudes, i.e., when going from t- to s-channel amplitudes all three helicities of the virtual 
photon contribute whether it is crossed or not. Interestingly, the dominant contribution for the 
kinematics we are interested arises if its helicity will flip, too. Hence, either in the s- or in the 
t-channel one can take that both photons have the same helicity in the center-of-mass frame. 
In the real world, one has to take more care of the rotation of helicity amplitudes, namely, a 
given s-channel helicity amplitude is expressed as a linear combination of t-channel ones, and for 
massive and/or virtual particles all helicities contribute. 

The partial wave expansion of t-channel helicity amplitudes is given in terms of Wigner d- 
matrices d^^O), which are essentially expressed by Jacobi polynomials [162]. Here the quantum 
number J refers to the total angular momentum, and fi (v) is given by the helicity difference of 
initial (final) state particles. In the DVCS process, y*h —> hj^*\ the twist-two CFFs, considered 
in this paper, appear in the helicity amplitudes in which both photons are transversely polarized 
and have the same helicity [24, 107]. As discussed above, since one photon is virtual, one can take 
that the photons in the t-channel process 7*7**) — > hh also have the same helicities. For the rest 
of this section, we take them to be —1, denoted as •y*^^ (the amplitude for 7*^ follows from 
reflection). Hence, for scalar (hh) and spin- 1/2 particles (e.g., h'h') the Wigner <i-matrices with 
/i = (photons) and v = {+1,0,— 1} appear in the partial wave expansion: 

Ki<-il ^ {hlhl} and d io for {^t}' ( 145 ) 

The SO (3) partial wave expansion of the t-channel helicity amplitudes yields those for the 
17 We are indebted to M. Diehl for critique and discussion on this subtle point. 
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crossed CFFs H® and H®, defined via Eqs. (5) and (6) by the replacement Pi — > —Pi and 
q 2 — > The hadronic tensor in the t-channel is obtained by the same substitutions from Eq. 
(3). The evaluation of the helicity amplitudes in the center-of-mass frame is straightforward. 
Taking into account the angular dependence arising from the spinor bilinears, 

.tf( fth „m»V(A) = ±^^, (1 46) 



(9i - *>) • (Pi - Pi) '""2 > " cosf) 

we can read off the 'effective' partial waves in the SO (3) expansion of the t-channel CFFs H® 
and Ti.® for massless hadrons: 

**>(cos0, 0, s®, Q 2 ) = f>J + l) 1±{ - 1)J fj(sV,#, Q 2 ) ,™ 6 2 A i(°) ■ ( 147 ) 

V 1 — cos^ V 



Here the upper (lower) sign refers to the parity even (odd) case, in which the total angular 
momentum is even (odd). The property of the Wigner matrices ensures then that the CFFs as 
scalar valued functions are invariant under reflection. The lowest partial wave that appears is for 
a total angular momentum J = 2 (1). Up to a phase factor, which is not indicated, the crossing 
relation for helicity amplitudes simply reads 

7i )± i/ 2; i i± i /2 ( S , t, Q 2 ) = TV^p^it®, s (t) , Q 2 ) , t {t) = a , s® = t (148) 

and leads with cos 9 = —l/rj to those for the CFFs: 

tf, A 2 , Q 2 ) = (cos = -1/rj, 0, s® = A 2 , Q 2 ) (149) 

where the photon virtuality asymmetry ■& — r)/£, cf. Eq. (9), is invariant under crossing. Combining 
this relation with the general form of the OPE (28), where the Wilson coefficients, as functions 
of are invariant under crossing, we can read off within our conventions the crossing relation for 
conformal moments [122, 98] 

Fjir), A 2 ) = T] j+1 Ff\cos9 = -l/rj, s® = A 2 ) . (150) 

Here F. are the conformal moments of generalized distribution amplitudes [24, 163]. Finally, a 
formal comparison of the crossed version of the partial wave expansion (147) with the OPE (28) 
allows us to identify the partial waves 

d J H {rf) oc -Jl==d J 0jl (6)\ coae ^ 1/ri , 2< J<j + 1 forj = {1,3,5, ■■•}, (151) 

vV - 1 

4(r?) oc -^^dlMUe^i/r,, l<J<j + l forj = {0,2,4,---}, (152) 

Vv 2 - 1 
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in the SO (3) expansion of conformal GPD moments Hj and Hj, where J is an even and odd 
number, respectively. The symmetry of the Wigner d- matrix ensures that these partial waves are 
always even polynomials in rj of order J — 2 or J — 1. Consequently, within odd (even) j we see 
that Hj {Hj) is an even polynomial in t], as it is required from time reversal invariance. The order 
of these polynomials is j — 1 and j for odd and even values of j, respectively. Here we omit further 
technical details and now discuss the partial wave amplitudes. 

5.2 Ansatz for SO(3) partial wave amplitudes 

First, we would like to illustrate with two examples the potential power of crossing, which allows 
us to consider the t-channel process rather than the s-channel one. For DVCS off a (pseudo) scalar 
particle, in which the photon helicity is conserved, we have one twist-two CFF Ti, which belongs 
to the parity even sector [121]. The expansion of the corresponding conformal GPD moments is 
done in terms of the Wigner <i-matrix 

V J+1 dio(0)\cose=-i/ v = v H1 PA-l/v) with < J < j + 1, and J, (j + 1) - even, (153) 

which is a polynomial in r\ of order j + 1 expressed in terms of a Legendre polynomial Pj (cos 9) = 

d, Q (9). The meaning of the quantum numbers in the crossed channel is obvious: two photons, 

e. g., travelling oppositely along the z axes, form a two-particle state with zero magnetic quantum 
number, projected on the z axes, and producing a pair of (pseudo) scalar hadrons, cf. Fig. 3(a). 
The total angular momentum of the initial state entirely transfer to the orbital angular momentum 
of the hadrons and so the latter is equal to J . The minimal value of the orbital angular momentum 
J , determined from the magnetic quantum number of the photons, is zero. Let us have a closer 
look at the partonic subprocess, which is thought of as the production of a quark- ant iquark pair 
state, labelled by the conformal spin j + 2, and travelling close to the light cone in different 
directions. Since at leading twist-two, only chirally even operators appear, helicity conservation 
holds at short-distances. Namely, the quark and antiquark have opposite helicities and, obviously, 
their spins point in the same direction, see Fig. 3(b). The conformal spin j + 2 must be larger 
than J+ 1 and is an odd number in the parity even sector. Hence, all odd conformal partial waves 
with J — 1 < j contribute to the SO (3) partial wave amplitude with the total orbital angular 
momentum J. However, the quark-antiquark pair is also bound by confinement and so they also 
represent an intermediate mesonic state with spin J, which consists of J — 1 units of the orbital 
angular momentum and one unit of the total angular momentum arising from the aligned spins 
of the quark and the antiquark. This mesonic resonance then decays into two scalar hadrons. 
The quark (antiquark) must be combined with an antiquark (quark), which is picked up from the 
vacuum. To not alter the quantum numbers of the vacuum, however, a quark- antiquark pair must 
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be picked up with opposite spin (magnetic quantum number), as shown in Fig. 3(b). Suppose that 
in one hadron the spins of the quark and antiquark are opposite, then in the other ones they are 
aligned. This would be in contradiction with a naive quark model picture and it can be resolved 
only if one of the produced quarks flips his helicity due to nonperturbative effects. 

Let us come back to DVCS off a nucleon, where again s-channel helicity conservation for 
the photons is required. Then the four CFFs TC, - ■ ■ ,£ can be expanded with respect to the 
t-channel partial waves d$ ±l and g?q . However, in our fictional world of massless particles, the 
helicity nonconserving quantities £ and E should vanish and so d^ is absent. Hence, we are 
left only with final states which have opposite helicities, i.e., 7i and 7i, which are expanded with 
respect to the partial waves (151,152). We can again employ our intuitive picture, and as in the 
intermediate resonance before, the quark and the antiquark form a spin-one state with the orbital 
angular momentum J — 1. By picking up two quarks and antiquarks from the vacuum the final 
spin-1/2 nucleon and antinucleon are formed and one would expect that they carry the helicity 
of the produced quark and antiquark, respectively. The corresponding conformal GPD moment 
for odd conformal spin j + 2 is now a polynomial of order j — 1. In the case that the helicity 
of the produced quark or antiquark is flipped owing to nonperturbative effects, the helicities of 
the produced nucleons will be the same. The partial wave g?q appears in the conformal GPD 
moments Ej and also Hj, which are polynomials of order j + 1. 

Let us emphasise that we have therefore observed a relation between the helicity flip and 
the order of conformal GPD moments. In general, the conformal GPD moment Hj is an even 
polynomial in rj, which is of order j + 1 for odd j. The minimal value of J is now given by 
the magnetic quantum number of the final state and so it is one. We remind ourselves that in 
the partial waves (151) a term of the order j + 1 in rj is absent and so it cannot appear in Hj. 
However, if we would allow for helicity nonconserved quantities, such a term arises, owing to the 
partial wave c?q , which will then appear in Ej and Hj, see the helicity representation for GPDs 
in Ref. [106]. We conclude again that the restoration of full polynomiality, namely, up to order 
(j + 1) for odd j, appears owing to the nonperturbative interaction in which the helicity of the 
parton is reversed. In other words, the breaking of chirality is encoded in the conformal GPD 
moments, which for given odd conformal spin j + 2 are polynomials of order j + 1. Certainly, the 
highest possible term in rj naturally arises as a part of the SO (3) partial waves. At this stage we 
see no reason to treat such terms in a special way and to collect them into a separate so-called 
D-term [120], which was introduced to complete the common spectral representation of GPDs. 
As mentioned in Sect. 3.2.2, we suggest to cure the spectral representation, e.g., as proposed in 
Ref. [121], rather than to add a D-term. 

The partial wave expansion we wrote down for GPDs was borrowed from the two-photon 
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fusion into two hadrons. In the case of hard vector meson production, the leading twist-two 
contributions arise from a longitudinally polarized photon and vector meson. It is easy to name 
the Wigner matrices, which are the same as for DVCS, namely, we have <io±i an d doo matrices 
in the SO (3) expansion of the conformal GPD moments Hj and Ej. The former ones, given by 
sin (6 , )Cj{ 2 1 (cos 9)/ a/ J( J + 1), are expressed by the Gegenbauer polynomials with index v = 3/2 
and the latter ones are the Legendre polynomials, or, if one likes, Gegenbauer polynomials with 
index v = 1/2. Again, the helicity nonconserved quantity Ej is formed from g?q , while Hj contains 
c/qj and an admixture of the cIq waves. If we now replace the vector meson by a pseudoscalar 
ones, the same partial waves appear for the parity odd quantities Hj and Ej. 

We realize that the SO (3) partial wave expansions of GPDs are universal, i.e., are the same 
for the considered processes. Such an expansion has several advantages, e.g., for the analytic 
continuation of even to odd j values and allows for a simple implementation of the normalization 
at j — 0. It also leads in a natural way to a term of order r]^ +1 in the polynomial Ej for odd 
values of j. Moreover, it is convenient to have for helicity conserved quantities the same functions 
that also appear in the conformal SO(2,l) representation. In a conformally invariant world there 
would appear only one SO (3) partial wave with J = j + 1. 

Inspired by the hadronic view on the t-channel scattering process, see Fig. 3(c), we now propose 
an ansatz for the partonic partial wave amplitudes that appear in the SO (3) expansion of the 
conformal GPD moments. Thereby, we rely on the Regge description of high-energy processes, 
which states that the high-energy behavior of the s-channel process "f^h — > h'j is dominated by 
linear Regge trajectories a(t) of pomeron and meson exchanges in the t-channel. The strength of 
the photon-photon-to-meson (pomeron) coupling is contained in a vertex factor fjj that depends 
on the conformal spin, too. The conformal spin j + 2 is considered as a variable conjugated to the 
partonic momentum fraction x. Thus, as in the case of mesonic distribution amplitudes that are 
expanded with respect to the Gegenbauer polynomials C^ 2 (x), fjj can be viewed in the partonic 
language as a probability amplitude for finding a quark-antiquark pair state with conformal spin 
j + 2 inside of a meson with given spin J. Furthermore, the partial wave amplitudes also contain 
the propagator l/(m 2 (J) — t) oc 1/(J — a{t)) of the exchanged particles, as well as the impact 
form factor, describing the interaction with the target. These form factors will be modelled by a 
p-pole ansatz, i.e, monopole {p = 1), dipole {p — 2), and so on, with a J-dependent cutoff mass 
squared M 2 (J). 

All together, we propose the following ansatz for the helicity nonflip conformal GPD moments 
in terms of partonic partial waves amplitudes, which we write down as a sum over the angular 
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momentum (we employ the fact that the polynomials are even): 

M(.|A',fl- V hi ^ [^"M for,-! * (154) 



l,odd 



where for odd (even) j the sum runs over even (odd) J. Here <^o iC^) ^ 7 ? J ^o / v) / \/iv 2 ~ 1) 
is the 'crossed version' (151,152) of the Wigner matrix, where the rotation matrix of the spinor 
bilinears is taken off. It is simply our partonic toy model (49), defined for complex valued j in 
Eq. (62): 



dti j+ \v) = ; j w & i j,j 1 - ^— j . (155) 



r(3/2)r(3 + j) r f-jj + 3 



2 1 +^T(3/2 + j) ' z 1 V 2 



n 



2'/ 



For our later convenience, the crossed rf-matrices are normalized in the forward limit to one: 

limdi„(i7)=l. (156) 

The conformal GPD moments (154) are even polynomials in 77 of order j — 1 or j, as required. 
We also note that they are build from even polynomials ^^(cos^) oc cos 6Cf^ 2 (cos 6), where 
/ = {0, 2, • • ■ , j — 1} or {0, 2, • ■ ■ , j}, with eigenvalue +1 under parity transformation. For helicity 
nonconserved quantities Ej and Ej , analogous ansatze can be written down in terms of Legendre 
functions. Thereby Hj will get an admixture from d^ partial waves, too. 

5.3 Modelling of conformal GPD moments 

It is beyond the scope of this paper to present a thorough study of realistic ansatze for GPD 
moments. We would rather like to convince the reader that the proposed parameterization gener- 
ically works and then use some simplified version for our numerical studies. One important aspect 
is the skewness dependence and its approximations. Another, new one, is the implementation of 
lattice results. 

The problem of how different approximations of the skewness dependence will affect the size 
of the corresponding CFF will be investigated within a toy model. It is similar to Eq. (49), which 
describes the Compton scattering process at tree level or, in other words, within a noninteracting 
parton picture. To make it somewhat more realistic, we multiply it with the generically valid 
Mellin moments of an unpolarized valencelike parton density 35a; _1 / 2 (l — a;) 3 /32 (here normalized 
to one): 

t _ r(i/2 + j)r(9/2) £ +1 _ 
M * Uj "r(9/2 + J )r(i/2)^ (i5 ° 
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Figure 4: The imaginary (left) and real (right) part of the CFF 7i(£, A 2 = 0) arising from the toy 
ansatz (157) for conformal GPD moments within different approximations (159): exact (solid), 
next-to- leading (dash-dotted), leading (dashed) SO (3) partial waves and £ = (dotted) term. 

Here the partial wave d^irj), given in Eq. (155), is normalized to one for vanishing skewness rj. 
For the sake of illustration, we now expand Hj° 7 (^) with respect to Legendre polynomials: 

d j+1 (n) - $ M + V -T 1 + ( V -T l ~ 2 r(/ + 3/2TO + 1) j' (n) (158) 

d H W ~ SoW + 4 2^ 2 [ 2 ) r{j + 3/2)r(/ + qdaflW > U58j 

where d l 0fi (rj) = (rj/^YV^^^ + l ) p i( l /v)/r{l + 1/2). The first term on the r.h.s. is the leading 
partial wave, while all other partial waves, contained in the remaining sum, are suppressed by 
powers of rf. We will now study the numerical deviation in the CFF that is induced by an 
approximation of the partial wave d 3 ^ (rj). For this purpose we make three choices: we drop the 
skewness dependence altogether, take the leading partial wave, and include also the next-to-leading 
ones: 

4 + -o = 1 , 4 + -lo = %,M > 4 + -nlo = %,M + T^-^^Co 2 ^) • (159) 

Note that some care is needed in the truncation of exact partial waves. The approximated partial 
wave cij^NLO differs for j = 1 (and also j = 0) from the exact one by rj /27 (3/r/). This is induced 
by a pole in ^ooC 7 ?) (^ool^)) ^ na ^ cancels the zero at j = 1 (or j = 0) in the expansion coefficient. 
Finally, this leads to an addenda in the real part of the CFFs that must be subtracted. In our toy 
example the subtraction term for the next-to- leading approximation reads —5/27r]. 

The outcome for the imaginary and real part of the CFF 7i* oy (£) is displayed in the left and 
right panel of Fig. 4, respectively. In the next-to-leading approximation (dash-dotted) one can 
hardly see a difference to the exact CFF (solid). If we take only the leading SO (3) partial wave 
(dashed), we realize that the deviation from the exact CFF (solid) is small over the whole £ 
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region, in particular for the imaginary part. Furthermore, making the expansion in powers of £ 2 
and retaining only the leading £ = term (dotted) yields the imaginary part that approaches the 
exact result already for £ < 0.7. However, it has an unrealistic feature that it does not vanish 
in the limit £ — > 1, as it should. Concerning the real part, both of these approximations start to 
digress from the exact result for £ > 0.3, where the deviation is larger for the expansion in £. 

To summarize, the expansion with respect to the angular momentum looks promising and 
works rather well in the case that the leading pole factorizes. We have also found in the expansion 
with respect to Legendre polynomials, which was used in Ref. [91], that the minimal version of 
the dual model, i.e., taking leading and next-to-leading partial waves, yields also an astonishing 
agreement with the exact result. We have also observed that in the experimentally accessible 
kinematical region the expansion in £ practically coincides with the exact result, in particular 
for the imaginary part. Nevertheless, there could be a drawback if the (leading) Regge poles 
are nonfactorizable as it is the case in our model ansatz (154). The next-to- leading term in the 
partial wave expansion generates an 'artificial' pole, e.g., l/(j — 1 — a(A 2 )), that is situated on 
the r.h.s. of the leading Regge pole l/(j + 1 — a(A 2 )). Such a pole leads to an addendum in the 
Mellin-Barnes representation for CFFs with the same small £-behavior as the leading pole. Thus, 
we expect that the normalization of the CFFs, e.g., for A 2 = 0, is governed by all terms in the 
partial wave expansion. A closer look at this potential problem should be given somewhere else. 
For the time being, we rely on the leading term in the expansion with respect to rj, 

j + 1 - «(A 2 ) (1 - j^)p 



Fjir,, A 2 , /i 2 ) = - ; ; J+1 ^ A2W1 Aa ,„ + 0( V 2 ) , (160) 



which we will use in our numerical studies. 

Next we fix the normalization of conformal GPD moments at A = 0: 

Fj(v = 0, A 2 = 0,/i 2 ) = - Jp±< . (I6l) 

I + 3 ~ a(0; 

The conformal moments of helicity nonflip GPDs Hj and Hj are then reduced to the Mellin 
moments of unpolarized (q) and polarized (Aq) parton densities, 

: = f dxx>q(x lf M>) and kj+1 ' 3 ' = [ dxx^Aq(x, Mo ) , (162) 

1+3 -a(Q) Jo 1+3 -a{0) J 

respectively. They are parameterized with the guidance of Regge phenomenology, determining its 
small x behavior to be x~ a ^ , whereas counting rules suggest their large x behavior, parameterized 
as (1 — x) 13 for x — > 1. Corresponding to Eq. (162), the generic ansatz for parton densities in the 
x space yields the following Mellin moments, e.g., 

g(x,/i 2 ) = Nx- ao (l-xf => h j+1J = N{l+j-a{0))B{l-ao+j,l3 + l), (163) 
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where B(a,b) = T(a)T(b)/T(a + b) is the Euler beta function. Note that the intercept «o of the 
Regge trajectory occurring in our ansatz and that obtained from a given fit of parton densities 
must agree. Hence, the leading pole in the Mellin moments (163) is cancelled and replaced by the 
corresponding trajectory 18 . Certainly, a more realistic ansatz for parton densities, to be used in 
global fits, could be obtained by a linear combination of such building blocks. In practice, either 
one can take the Mellin moments of one of the standard parameterizations of parton densities 
[3, 4, 5, 6, 7, 8, 9, 10] or, if they are plagued by larger errors or theoretical uncertainties, one can 
perform a simultaneous fit of exclusive (e.g., DVCS) and inclusive data. Here one has to bear in 
mind that parton densities are scheme- dependent quantities. 

Further constraints for the conformal moments arise at j — 0. There they are given by the 
so-called partonic form factor 

^ofaA-) = 1 _^ A 1 (164) 



After adjusting the flavor quantum numbers, the partonic form factors coincide with the measured 
elastic form factors of the proton. The free parameters left, i.e., the cut-off mass Mq and the power 
behavior p at large A 2 can be taken from a fit to experimental data, where a refinement of our 
parameterization might be necessary. We note that form factors are scheme independent, except 
the axial one in the flavor singlet sector. Another advantage of the parameterization (154), which 
we will not use here, is that the partonic partial wave amplitudes fjj are related to physical ones, 
denoted as fj(W 2 , Q 2 ), in the t-channel. The physical amplitudes fj are given as a series of 
partonic ones, where the sum runs over the conformal spin j + 2. Note that such a relation depends 
on our scheme conventions, too. 

Finally, plugging the normalization (163) into the ansatz (160), we end up with the following 
simplified GPD parameterization: 

Fjir,, A 2 ,/i 2 ) =NB(1 - a(0)+j,P+ l) j±±^± + 0( V 2 ) . (165) 

We remark that for large |A| 2 , the counting rules predict a power-like falloff of form factors as 
(l/|A| 2 ) ns , where n s is the number of spectators, while the large x behavior of parton densities is 
(l-x) 2n °-\ i.e., 

p = n s - 1 , (3 = 2n s -l. (166) 

An inclusive-exclusive relation between the unpolarized DIS structure function W 2 an d the elec- 
tromagnetic form factor F 2 has also been derived by Drell and Yan within a field-theoretical model 



18 Remaining nonlcading poles at j = {—2 + uq, —3 + cto, ■ • ■ } arc considered an artifact of the parameterization 
and appear as subleading contributions in the CFFs. 
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Figure 5: The A 2 dependence of the conformal GPD moments (167) for valence quarks (left) 
and the spin dependence of the effective dipole mass (170) (right). In the left panel the solid line 
displays the dipole fit (169) to the experimental data and the dash-dotted one is our ansatz (167) 
with j = 0. The dashed and dotted lines show the changes with j = {1, 2, 3} for our ansatz and 
a dipole fit within the dipole masses (170), where AM 2 = (M p ) 2 . In the right panel we show the 
effective dipole masses (170) as a function of spin within AM 2 : (dashed), M 2 (dash-dotted), 
(2M P ) 2 (solid). The dotted line has a slope that arises from AM 2 = 2(2M P ) 2 and is compatible 
with lattice measurements [165] for the flavor nonsinglet combination u — d in the heavy-pion 
world, which was linearly extrapolated to the physical pion mass. 

that accounts for the dynamics of partons [164]. This result coincides with the counting rules for 
valence quarks. However, one should be aware that for sea-quarks and gluons, these counting rules 
might be modified; for a discussion see, e.g., Ref. [2]. For simplicity, we will not account for that 
here. 

Finally, we would like to demonstrate that our ansatz (165) can be easily adjusted to the 
experimental data on electromagnetic form factors and, moreover, that it is well suited to include 
lattice data. In particular, present lattice measurements [35, 36, 37, 38, 39, 40] give insight into 
the functional change of the A 2 dependence with increasing conformal spin. This dependence 
arises from two sources: the Regge trajectory and the impact form factor. As an example, we 
consider the conformal GPD moments of valence quarks, where the Regge trajectory is generically 
correctly described by a(t) = a(0) + a't with the intercept a(0) = 1/2 and the slope a' = 1 GeV -2 . 
For two spectators the counting rules state that the impact form factor is a monopole, i.e., p = 1. 
For its cutoff mass we naturally choose two times the proton mass M Q ~ 2M p = 1.88 GeV. Hence, 
our conformal GPD moments are written as 

Hf( V , A 2 , /i 2 ) cc \ 2 \ 2 + 0( V 2 ) . (167) 

J 3 
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Here we have introduced the monopole mass squared m 2 = (1 — a(0) + j)/o! = (1/2 + j) GeV 2 , 
which arises from the mesonic Regge trajectory. We might also assume that the spin dependence 
of the monopole mass squared M 2 is linear 

Mf = M 2 + AM 2 j . (168) 

The generic ansatz (167) leads to a satisfying description of the electromagnetic proton form 
factor. In Fig. 5 (a) we confront this ansatz (dash-dotted) with the dipole fit (solid) of the 
electromagnetic proton form factor F\, 

Fl (Q 2 = -A 2 ) = / ^^f 606 ^, , m 2 dipol = 0.71 GeV 2 , (169) 



and realize that they fairly agree. We also display the A 2 dependence of the conformal moments 
within the choice AM 2 = (M p ) 2 . A larger value of AM 2 is compatible with the slope measured 
on lattice in the heavy-pion world, as shown in Fig. 5 (b) (dotted line). However, the intercept 
differs (in fact, the lowest moment does not describe the Fi(Q 2 ) data). The masses are extracted 
from a dipole fit to the lattice data. To compare with our ansatz, we calculated the 'effective' 
masses 

= 2 ( i +i - n( o) + wvxm) 1 • (170) 

appearing in a dipole fit. For |A 2 | < 1 GeV 2 such a 'refitting procedure' only weakly modifies the 
GPD moments, compare dashed and dotted lines in Fig. 5 (a), and so this procedure is justified 
to some extent. 



6 Perturbative corrections of the DVCS cross section 

This section is devoted to the numerical analysis of radiative corrections to CFFs. We will con- 
centrate on the CFF H, since it is the dominant contribution in most of the DVCS observables. 
The three remaining twist-two CFFs are usually suppressed by kinematical factors [68] . To reveal 
these CFFs from experimental data, it is therefore crucial to understand the theoretical uncer- 
tainties of 7i . Our findings can be qualitatively adapted for the helicity nonconserved CFF S, 
which enters the parity even sector, too. In this sector we face a peculiarity that is related to 
the appearance of the pomeron trajectory in high-energy scattering. Technically, it shows up as 
an essential singularity of the evolution operator at j = 0, see anomalous dimensions (108) and 
(109). Such a singularity is absent from the parity odd sector. Hence, our results in the small £ 
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region are not directly applicable to the study of radiative corrections of the two remaining parity 
odd CFFs H and £. 

Based on our model (154), we introduce in the next section a simplified generic ansatz for 
conformal GPD moments, which will serve us in our numerical studies of the experimentally ac- 
cessible kinematical regions. In Sect. 6.2 we shortly discuss the features of flavor singlet and 
nonsinglet CFFs to LO accuracy. In the following two sections we then elaborately analyze the 
size of radiative corrections to NLO and, finally, to NNLO accuracy. Thereby, independently of 
the considered order and scheme, we take the same conformal GPD moments which thus leads to 
different CFFs. However, note that CFFs are physical observables and do not depend on our con- 
ventions. Therefore, conversely, the conformal GPD moments revealed from a measurement of the 
physical CFFs will depend on both our scheme conventions and the approximation. Nevertheless, 
the results of our analysis give us a measure for both the reparameterization of the GPD ansatz 
needed to compensate convention change and for the convergency of the perturbation theory. 



6.1 A simplified generic ansatz for conformal GPD moments 

The kinematics of interest can be restricted to rj = £ < 0.5, i.e., the Bjorken scaling variable 
XBj = 2£/(l + £)is bounded by 2/3. For simplicity, we do not resum the partial waves and rely 
on the leading term in the t] expansion (165). 

Let us specify the other parameters in the ansatz (165) for a valence-like helicity non-flip GPD 
H. The leading meson Regge trajectory is generically given by a(t) = a(0) + alt with a(0) = 1/2 
and a' = 1 GeV -2 . For two spectators the counting rules (166) state that p = 1 and (3 — 3. For 
the cut-off mass of the impact form factor we choose (168) with M = 2M P = 2AM = 1.88 GeV. 
Hence, we have for our GPD moments: 

A 2 , A« = B pg^\ V rAr + , ( m ) 

^W^j i_ (l+2i)GeV2 l ~M^+D 

which is normalized to one for j — and A 2 = 0. 

The conformal GPD moments with the flavor nonsinglet combination (A. 7) or (A. 8), which 
is also relevant for DVCS, is built from valence and sea quarks. For four active quarks and the 
SU(2) flavor symmetric sea, the sea quark part arises from the difference of charm and strange 
( ant i-) quarks, cf. Eq. (A. 8). Suppose we are at the charm threshold and the charm sea is generated 
dynamically. Then, essentially, only the (anti-) strange quark counts. We might assume that the 
u, d, s antiquarks have the same conformal GPD moments and so the breaking of SU(3) flavor 
symmetry is described by one single parameter R s /u, defined as the ratio of s to u antiquarks. For 
the purpose of illustration we consider two alternative cases: one without and one with sea quark 
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admixture: 



NS ^,A a ,^)= J H7 a (^A 2 > ^), (172) 
™ Hj (rj, A 2 , = Hf( v , A 2 , /, 2 ) - J^-iq-iTl, A 2 , ^) • (173) 

For the valence quark content we rely on Eq. (171) and the sea quark GPD moments are specified 
below in the ansatz (175) with iV sea = 4/15. For the SU(3) flavor symmetry breaking parameter 
we choose the often used value Rg/u = 1/2. 

In the flavor singlet sector we need an ansatz for both quark and gluon conformal GPD mo- 
ments: 

Hj(r}, A 2 , /i 2 ) = ( H J\ (77, A 2 , vl) , Hf = Hp + Hp* + Hf* = H s ™ + 3Hf . (174) 

The singlet quark combination consists of sea and valence quarks, cf. (A. 7) and (A. 8), and we 
employ isospin symmetry to express the latter within the ansatz (171). Moreover, we rely on the 
counting rules (166), i.e, /3q = 5,^g — 2,/3 sea = 7,p sea = 3, and take the same cut-off mass for the 
impact form factors as for the valence quarks before. Hence, analogously to Eq. (171), the generic 
ansatz (165) then reads: 

^ - "- Vt KV^, ri * y + ° ( " 2) ' (175) 

ttG( A 2 2n m B ( l ~ «g(0) +J,6) 1 1 2 

where = M| ea = M 2 (4 + j) and m 2 = (1 — a(0) + j)/a' is expressed in terms of the intercept 
and slope of the Regge trajectories, specified below. 

The leading trajectory arises now from the 'pomeron' exchange. We remind that in deeply 
inelastic scattering the structure function F 2 ~ (1/xbj) a ^ 2 ^ grows with increasing Q 2 . Here the 
exponent is governed by the intercept of the Regge trajectory A = a(0) — 1, which is, in the 
language of Regge phenomenology, that of the soft pomeron (for Q 2 — > 0): 

ap(t) = a P (0) + apt a P (0)~l, a' F = 0.25 . (177) 

However, in hard processes the trajectory will effectively change due to evolution, which differently 
effects the behavior of quark and gluon parton densities. In particular, the value of a(0) increases, 
while that of a' decreases with growing resolution scale Q 2 , e.g., see Ref. [166]. In the flavor 
singlet sector the size of radiative corrections and the strength of evolution crucially depends on 
the effective pomeron parameters, see for instance the variation of NLO corrections obtained in 
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[85]. In our numerical studies we shall consider two scenarios in which the radiative corrections 
to NLO accuracy are respectively small and large. It is known that the corrections to the quark 
sector are rather stable, while the main uncertainty arises from the gluons, which enter in the 
perturbative description of CFFs at NLO. Small and large NLO corrections can be obtained by 
choosing a softer and harder gluon, respectively: 

"soft" gluon: N G = 0.3, a G (0) = a sca (0) - 0.2 , (178) 

"hard" gluon: N G = 0.4, a G (0) = a sca (0) + 0.05 . (179) 

Furthermore, we choose a realistic value for a sca (0) = 1.1 and a' sea = a' G = 0.15, all at the input 
scale Ql = 2.5 GeV 2 . 

We remark that we normalize the sea quark (175) and gluon (176) moments at j — 1, so A^ ea 
and Nq give the amount of momentum fraction carried by the considered parton species. Because 
of the momentum sum rule (A. 20), valid in the forward kinematics, we have the constraint 

N G + iVsca + / dx x [uy + dy] (x) = l. (180) 
Jo 

Within our toy ansatz for valence quark moments (171) we have for their momentum the generic 
value 1/3 and so iV sca = 2/3 — Nq. We note, however, that the separate contributions will change 
during the evolution. In the asymptotic limit Q — > oo, the evolution equation tells us that 
Nq = 4CV/(4Cj? + n/), i.e., that more than 50% of the longitudinal proton momentum is carried 
by gluons. As it is experimentally verified, at a scale of a few GeV 2 the gluons already carry about 
40 % of the momentum. 

6.2 Compton form factors to LO accuracy 

The CFFs ty-L and NS 7i are evaluated from the specified conformal GPD moments by the Mellin- 
Barnes integral (58). Essentially, we employ here the same technique that is well established in 
deeply inelastic scattering, only the integrand is now more intricate. The integral does not depend 
on the integration path, going from c — zoo to c + zoo, as long as we do not cross any singularities, 
see Fig. 2 (a). In practice we use this property to get closer to the leading singularity, lying left of 
the integration path. Since the conformal GPD moments rapidly decrease with growing conformal 
spin j, we can in practice cut the integration path to a finite length. Moreover, the convergency 
can be improved by a separate rotation of the integration path in the upper and lower half-plane 
in such a way that along the path the real part decreases, as shown in Fig. 2 (b). Nevertheless, 
it is always a good idea to exercise proper care in the choice of the integration path and to check 
the numerical accuracy. Once this is done for given conformal GPD moments, the numerical 
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treatment is simple, fast, and stable even in NNLO. We use two different codes for the numerical 
evaluation, one specifically written in FORTRAN and, alternatively, the integration routine from 
MATHEMATICA. 

For the evaluation of the CFFs to LO accuracy, we use the flavor nonsinglet Wilson coefficients 
(91a) and (92) as well as the singlet ones (101a) and (102). The Q 2 evolution is governed by the 
flavor nonsinglet operator (116) and singlet one (121)— (123), approximated to LO. We equate the 
factorization scale with the photon virtuality: fi 2 = Q 2 . The various ansatze, given in Eqs. (172), 
(173), (175), and (176), are taken at the input scale <2q = 2.5 GeV 2 . The running coupling in LO 
approximation is normalized to cx s (Qq = 2.5 GeV 2 ) jix = 0.1, where the number of active quarks 
is four. 

We found it useful to describe the complex valued CFFs in polar coordinates, 

H(£, A 2 , Q 2 ) = | A 2 , Q 2 ) | exp {zp(£, A 2 , Q 2 )} , -ix < <p = arg (H) < tt , (181) 

rather than Cartesian ones. This avoids a discussion of radiative corrections in the vicinity of 
zeros, appearing, e.g., in the real part of certain CFFs. Moreover, the polar coordinates reveal a 
simple shape of CFFs in their functional dependence on £ and Q 2 . This is demonstrated in Fig. 6, 
where we display the modulus (left) and the phase (right) of the pure valence nonsinglet (up) and 
the "hard gluon" singlet (down) CFF, respectively. For the former one we choose the kinematical 
region that covers the phase space of present fixed target experiments, namely 0.05 < £ < 0.5 
[0-1 < ^Bj < 0.65] and 1 GeV 2 < Q 2 < 10 GeV 2 . For the latter one we also include the phase 
space that is explored in collider experiments, i.e., 10^ 6 < £ < 0.5 [2 ■ 10~ 6 < xej < 0.65] and 
1 GeV 2 < Q 2 < 100 GeV 2 . The simple shape of both the modulus and the phase is perhaps a 
more general feature of CFFs and should not be considered as an artifact of our approximation 
[85]. It is for 1 GeV 2 < Q 2 < 10 GeV 2 almost independent on the photon virtuality. In Sect. 6.3.1 
we will have a closer look to the evolution. Here we remark only that both the moduli and phases 
of CFFs are rather "planar" in fixed target kinematics. With increasing l/£ the evolution starts 
strongly to affect the modulus in the singlet sector. It grows with increasing scale Q 2 , namely, 
in such a way that the resumed logarithmical scaling violations lead to a power like change of 
its l/£- dependence. This feature is well known from unpolarized DIS and in agreement with 
experimental DVCS results [65, 66, 67]. The phase is much less affected by evolution and in the 
small £ region it is nearly independent of £. It monotonously decreases from its input value to be 
changed by less than 25% at Q 2 = 100 GeV 2 . 
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Figure 6: The scaled moduli (a,c) and phases (b,d) of the flavor nonsinglet (up) and singlet 
(down) CFF Ji are displayed to LO accuracy versus Q 2 and £ for fixed A 2 = —0.25 GeV 2 . At the 
input scale Qq = 2.5 GeV 2 we use the ansatze (172) and (179) for conformal GPD moments. The 
running coupling is normalized to a s (Ql)/n = 0.1, where the number of active quarks is set to 
four. 



6.3 Size of NLO radiative corrections: CS versus MS scheme 

Now we explore the radiative corrections to NLO accuracy and in particular the differences between 
the CS and MS schemes, separately for the flavor nonsinglet and singlet sector. Before we do so let 
us explain the relation between these schemes and define the quantities that serve us as measure 
for the size of perturbative corrections. 

The perturbative expansion of a CFF in the CS scheme, e.g., for the flavor nonsinglet case, 
might be structurally written to NLO as 



c (o) + ^i( C (i)cs + C (°)AW)\ £ {0) ® F cs + 0{a 2 s ) . (182) 



Here £^ denotes the evolution operator (116) in LO approximation and is the NLO cor- 
rection. The convolution symbol ® indicates the integration over the complex valued conformal 
spin j, which is for shortness not particularized, where the measure includes the appropriate nor- 
malization and £ dependence, see (58) and (91a). By construction, in the forward kinematics 
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this scheme is identical with the MS one, see normalization condition (75). However, for DVCS 
kinematics the form of Wilson coefficients and evolution operator in the MS scheme are already at 
NLO modified by off-diagonal, i.e., rj proportional, terms, see Eqs. (72) and (73). Because of the 
particularity of the DVCS process, where rj — £, both Wilson coefficients and evolution operator 
are finally modified by an infinite sum of terms, which appears then in front of the conformal 
GPD moments, cf. Eq. (141): 



C<0) + ^(C(1)MS + ^(0)^(1) + C (0) £(0) g, F MS + (lg3) 



Here © i^ 1 -* contains off-diagonal part of the evolution operator, where © symbolizes the 
summation over the conformal spin, cf. Eq. (142). Since the physical observable T is independent 
of our conventions, the conformal moments in both schemes are related to each other by a scheme 
transformation. At the input scale Q 2 = <2q, where = 1, = 0, and = 0, we might 
express this transformation by a finite factorization (or renormalization) constant 

F m =F m +^zU(r,)®F™, C^ = C^ m + C^®z^. (184) 

27T 

For (positive) integer conformal spin the triangular matrix z^{rj) = rf~ k z^ contains only off- 
diagonal entries, i.e., j > k + 2. Note that the change of Wilson coefficients is of course r]- 
independent, while the skewness dependence of the conformal GPD moments is altered, which is 
at least suppressed by a factor rj 2 . In particular, the two lowest GPD moments are untouched by 
the scheme transformation, while the (positive non-vanishing integer) moments are modified by 
a s and rj 2 suppressed contributions: 

F ms = F cs, F MS = F p 5 and F p = F cs + %LO(ri 2 ) for j = 3,4,---. (185) 

Strictly spoken, the truncation of the perturbative expansion also induces a discrepancy in the 
CFFs between these two schemes, which is beyond the approximation we are dealing with, i.e., of 
order a 2 . However, the conformal GPD moments revealed from a given data set will differ to order 
a s . In the following we study the NLO corrections in both schemes within the same ansatz. The 
resulting deviation in the CFFs can be viewed as a measure for the needed reparameterization of 
conformal GPD moments by altering their r, dependencies. 

We introduce now the quantities that we utilize as a measures of the scheme dependence and, 
foremostly, as indicators of the convergency of the perturbation series. It is natural to employ for 
this purpose the ratio of the CFF at order N p LO to the one at order N p_1 LO, where P = {0, 1, 2} 
stands for LO, NLO, and NNLO order, respectively: 



r (£, A 2 , Q 2 \Q 2 ) = K A , Q^|Qo)exp{i(TV(f, A , Q 2 \Q 2 )} . (186) 
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The phase difference 

^ = ar g(^LoJ ~ K) ( 187 ) 

is formally of order (a s /27r) p . If convergency holds, it diminishes in higher orders. Under this 
circumstance the ratio of moduli 



K p 



h nPlo 



■^N^-iLO v 1 

approaches one, i.e., its deviation from one vanishes, too: 

5 P K = K p - 1 ~ O (of) . (189) 

If we suppose that the perturbative expansion is not ill-behaved, then the radiative corrections 
should not overshoot the size itself of the CFF in a given order, i.e., 



with r < 1 . (190) 



If this inequality holds true, the phase difference (187) is geometrically constrained by the value 
of r, otherwise it is independent. More precisely, we have the upper bound 

|sin<5 p y?| < r \$ P( p\ < ^ r f° r r < 1 . (191) 

Moreover, the triangle inequality, applied to Eq. (190), constrains the variation of the modulus, 
namely, r > 5 p iT|. Employing the cos-theorem, we can appraise the radiative corrections (190), 
quantified by the ratio r, from the variation of the phase and modulus. It will turn out that in 
our analysis the phase difference \S p <f \ is always small, i.e., \5 p ip\ tx/2. Hence, we can rely on 
the expanded version of this theorem, 

r 2 ^ (5 P K) 2 + (1 + 5 P K)(5 P V ) 2 , (192) 

which gives us a simple form of the constraint among the variations of the modulus and phase as 
well as the size of radiative corrections. Note that in the case of a (very) small phase change the 
variation of the modulus is roughly estimated to be |5 p i^| ~ r. 

6.3.1 Flavor nonsinglet sector 

The flavor nonsinglet CFF NS 7i is straightforwardly evaluated from the two ansatze (172) and 
(173) by means of the Mellin-Barnes integral (58). The Wilson coefficients, needed for our NLO 
(P = 1) analysis, are listed for the CS and MS scheme in Sect. 4.2.1 and 4.3, respectively, while 
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Figure 7: The relative NLO radiative corrections (189) and (187) to NS 7Y for A 2 = (left) and 
A 2 = — 1 GeV 2 (right) are plotted versus £ for the moduli (up) and phases (down) within the CS 
(dashed) and MS (dash-dotted) scheme. Thin and thick lines refer to a valence-like ansatz (172) 
and one with a sea-quark admixture (173), given at the input scale Qq = 2.5 GeV 2 . We equated 
the scales fif — fi r = Q and take the normalization condition a s {Q^)/-n = 0.1. 

the relevant expansion of the evolution operator can be read off from Sects. 4.2.3 and 4.3. As said 
above, we combine the perturbative expansion of the Wilson coefficients with that of the evolution 
operator in a consistent manner, where the leading logs are resummed, cf. Eqs. (182) and (183). As 
long as it is not stated otherwise we equate the factorization fif = /i and renormalization fi r scales 
with the photon virtuality Q. As input scale for the conformal GPD moments we use as before 
Qq = 2.5 GeV 2 . This input scale serves us also to normalize the coupling constant a s (Ql)/n = 0.1, 
where its running is described by the exact numerical solution of the NLO renormalization group 
equation. 

In Fig. 7 we display the relative NLO corrections (189) and (187) in the CS (dashed) and 
MS (dash-dotted) scheme to the moduli (up) and phases (down) for the phase space of present 
fixed target experiments. Thereby, we employ the pure valence-like ansatz (thin lines), i.e., Eq. 
(172), and the one with a sea quark admixture (thick lines), given in Eq. (173). For the mo- 
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mentum transfer squared we take two extreme values, namely, A 2 = (left) and A 2 = — 1 GeV 2 
(right). Comparing the resulting radiative corrections for both choices, one realizes that the 
A 2 -dependencies only slightly influence their size. 

Due to the radiative corrections the phases (lower panels) increase by a small amount, except 
for the valence-like ansatz (thin lines) where we observe a decrease for smaller values of £. In any 
case, the absolute value of the phase differences \$ P (p\ does not exceed 0.17T rad. The influence of 
the NLO corrections on the moduli (upper panels) is more pronounced. Generally, they moderately 
reduce the moduli, however, for the sea quark admixture ansatz in the CS (thick dashed) there 
is a small increase at larger values of £. The modulus is more affected for the valence-like ansatz 
(thin) than for the one with a sea quark admixture (thick), while in the case of the phase difference 
the situation is reversed. 

Comparing the corrections in the CS (dashed) and MS (dash-dotted) scheme, one realizes 
that in the former scheme they are smaller for the moduli, while the phase differences are almost 
independent of the specific choice. This is in agreement with the findings of Ref. [101], where 
a slightly different ansatz has been chosen. As explained above, the difference between the two 
schemes originates from the skewness dependence. For positive integer conformal spin, we would 
count them in the conformal GPD moments as if effects, suppressed by a s /2ir,. However, as 
we also spelled out, in DVCS kinematics we should not use r\ = £ as an expansion parameter, 
since the change of Wilson coefficients is 77- independent. Indeed, roughly spoken, the moduli 
differences in both schemes are of the same order as the radiative corrections themselves and 
nearly £ independent. This observation should be understood as a warning that rf suppressed 
terms in conformal GPD moments perhaps cannot be simply neglected by formal rj 2 counting. 

We consider now the evolution to LO and NLO accuracy, where we compare the modulus and 
phase at a given scale with those at the input scale, quantified by the ratios: 

^■^-ffl- 1 ^*'® = ^{m)- (193) 

Note that in contrast to the definitions (187) and (189), in which the variation is denoted by 
S, here we do not compare the CFFs in different order, but rather measure the strength of the 
evolution within a given order. The evolution in the MS scheme is consistently treated, i.e., the 
mixing of conformal GPD moments is taken into account. We remark, however, that this effect is 
tiny and can be safely neglected. For instance, for the quantity 



\ NS H(Q 2 ) 




NS H dia (Q 2 ) 




N 


%(Q 2 )| 



where the superscript 'dia' stands for neglecting the non-diagonal parts in the anomalous dimension 
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Figure 8: The modulus (up) and phase (down) changes of the evolved CFFs NS TC are plotted for 
fixed Q 2 = lGeV 2 (left), Q 2 = 10 GeV 2 (right), and A 2 = -0.25 GeV 2 versus f in LO (dotted) 
and NLO: CS (dashed) and MS (dash-dotted) scheme. The ansatze and scale setting prescriptions 
are the same as in Fig. 7. 

NS 7jfc, we find in a broad range of £ and Q a value on the level of few per mil. The phase differences 
are negligible, too. 

In Fig. 8 we show the evolution effects for the same ansatze as before, in LO (dotted) and NLO 
for both the CS (dashed) and MS (dash-dotted) schemes. In the left and right panels we plot the 
quantities (193) at a scale Q 2 = 1 GeV 2 and Q 2 = 10 GeV 2 , respectively. For simplicity, we do not 
perform any matching at the charm threshold. As already observed for the radiative corrections 
at the input scale, see lower panels in Fig. 7, the phase differences in right [left] panels in Fig. 8 
are again rather small and lie for forward [backward] evolution in the interval — 0.067T ■ ■ ■ — 0.027T 
[0.027T • • • 0.067r] rad. The signs tell us that the phases decrease during the evolution to a larger 
scale. Radiative corrections amplify this effect, where the differences between both schemes are 
again tiny. Comparing the upper left and right panels, we see that for the valence-like ansatz in 
LO (thin dotted line) the evolution mildly affects the moduli, i.e., of the order of about ±5% at 
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the edges of the phase space 19 explored in the present fixed target experiments. 

Furthermore, the modulus crosses for low and high values of Q 2 (upper left and right panels) 
the zero line at almost the same point £ ~ 0.12 to LO accuracy (thin dotted). This means that 
the CFF in the vicinity of this point is nearly scale independent. As in DIS, the evolution predicts 
a decreasing modulus for £ > £o and increasing one for £ < £o with growing Q 2 . The value of £o is 
a function of A 2 and depends on the input. In our example it is shifted by radiative corrections 
to the right (thin dashed and dash-dotted lines). This balance effect arises from the fact that 
forward evolution suppresses (enhances) the large (small) momentum fraction region. 

The radiative corrections lead to an amplification of the evolution effects. Note that the 
Wilson coefficients now also depend on Q 2 and the perturbative corrections are getting smaller 
with increasing Q 2 . Now the variation of the CFF modulus reaches for the valence-like ansatz 
(thin lines) the =pl0% [T^%] (for small £) to ±8% [±6%] (for large f ) level in the CS (dashed) [MS 
(dash-dotted)] scheme. For the ansatz with sea quark admixture (thick lines) there is no crossing 
point with the zero line anymore and so the 'balance-point' is shifted outside of the discussed 
kinematical region. Also here the strength of evolution grows by radiative corrections and varies 
in the range from ±6% (small £) to ±15% (large £). 

As we realized, in both ansatze the NLO corrections within the MS (dash-dotted lines) and the 
CS (dashed lines) change the LO prediction (dotted line). The differences, caused by the scheme 
dependence, mainly arises from the different Wilson coefficients to NLO, which are multiplied 
with the LO evolution operator. The contribution of the diagonal part of the NLO anomalous 
dimensions is small, i.e, about one percent. As we spelled out above, the influence of the non- 
diagonal part, appearing in the MS scheme in the anomalous dimensions, is negligible. 

We would like to briefly confront evolution effects with experimental measurements from the 
Hall A experiment at Jefferson LAB [75], where scaling was reported. Within the lever arm 
1.5 GeV 2 < Q 2 < 2.5 GeV 2 , we find for the valence-like ansatz that the scaling violation due to 
the LO evolution is small, namely, the modulus of the CFF varies by about 1.2% for xbj = 0.36 
(i.e., £ = 0.22) and A 2 = —0.25 GeV 2 , while a sea quark admixture can lead to a change of 
up to 4%. In NLO the variation for the former ansatz is —0.5% [—2%] for the CS [MS] scheme, 
while for the later one we find 5.7% [4.3%]. If one naturally assumes that for xb] = 0.36 the 
valence components dominate, one might conclude that the observed scaling in Ref. [75] indicates 
the smallness of higher twist contributions. However, in modelling of GPDs one usually realizes 
that the role of sea quarks in the CFFs is more pronounced than in DIS structure functions. 

19 For —A 2 ~ 1 GeV 2 this evolution effect might increase to become of the order of ±10% for smaller values of 
£. This is caused by the shift of the leading meson Regge pole to the left, approaching the j : = — 1 pole of the 
anomalous dimension, see Eq. (95). 
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variation, order/ £ 


0.05 


0.1 


0.25 


0.5 




3.7 [-6.9] 


0.7 [-8.5] 


-3.8 [-9.5] 


-8.0 [-10.9] 


// f , NLO (CS) 


0.5 [-0.8] 


0.7 [-0.7] 


0.8 [-0.3] 


-0.1 [-0.8] 


fi f , NLO (MS) 


0.3 [-0.6] 


0.7 [-0.3] 


1.0 [0.3] 


0.4 [0.2] 


// r , NLO (CS) 


5.5 [1.3] 


4.9 [0.6] 


2.9 [-0.4] 


0.3 [-1.6] 


fi r , NLO (MS) 


7.9 [3.9] 


7.6 [3.2] 


5.7 [2.2] 


2.9 [0.8] 



Table 1: Variation (194) in percent of the nonsinglet CFFs NS 7i induced by separate factorization 
and renormalization scale changes from Q 2 /2 ■ ■ -2Q 2 . Here we used the valence-like [with sea quark 
admixture] conformal GPD moments (172) [(173)] and set Q 2 = 4GeV 2 and A 2 = -0.25 GeV 2 . 



Interestingly, it has been argued that indeed the leading Regge trajectory essentially contributes 
to DVCS even in the valence region, in contrast to DIS [167]. We only like to point out here that 
even for fixed target kinematics a detailed view on scaling breaking effects for the net contribution 
to CFFs is necessary and that it might be used to constrain the GPD ansatz. 

So far we have considered only the scale setting prescription fif = \i r = Q. Change obtained 
by choosing another prescription within the same input scale Qo is often considered as an estimate 
for the possible size of higher order corrections. Let us have a closer look at this point of view. 
A change of \x T modifies the size of the Wilson coefficients by formally inducing a (3q proportional 
contribution a 2 /(27r) 2 that is multiplied with the NLO correction to the Wilson coefficients them- 
selves, see Eq. (91c). A modification of fif essentially corresponds to the difference between the 
expanded and non-expanded version of the evolution operator which is formally also of higher 
order in a s . However, besides non-leading log terms, it contains also leading ones, e.g., propor- 
tional to a 2 /(27r) 2 In 2 (//// Qo)- Certainly, whenever a new entry appears in the next order that is 
completely independent of these quantities then the rough higher order estimate can fail. 

Let us explore these estimates in more detail by employing the definitions 

x \H(Q 2 \fi 2 = 2Q 2 )\-\H(Q 2 \fi 2 =Q 2 /2)\ 

\H{Q 2 \ti = Q 2 )\ ' 1 ' 

where fi\ is the factorization (i = f) [renormalization (i = r)] scale and the renormalization 
[factorization] scale is fixed to be Q 2 . At LO we can only change the prescription for the ambiguous 
factorization scale setting. As one can read off from Tab. 1, the scale uncertainty to LO goes from 
about 4% [—7%] to —8% [—11%] with increasing £ for the valence-like [with sea quark admixture] 
ansatz. Comparing with the upper right panel in Fig. 8, we realize that this uncertainty reflects 
nothing else but the LO evolution itself. As we have expected, these numbers are not correlated 
to the perturbative corrections. At NLO the factorization scale dependence is drastically reduced 



64 



-80 
0.2 

a 

J 



~ i — i i i 1 1 mi 1 — i i i 1 1 mi 1 — i i i 1 1 mi 1 — i i i 1 1 mi 1 — i r 



— -20 h 

S - 40 ------ 

^ -60 



^ 



A 2 = 



J I I I I I I I I l_L 



-0.1 



J 1 


1 1 1 


lllj 1 1 1 1 1 INI 1 II 


l l l ll| l ll 


1 1 1 ll| 1 II 






NLO, CS 


A 2 - 


= : 






■ — ■ NLO, MS 










y 


y ' : 
















7*<* 




1 I 


1 1 1 


"""I 1 i " " 


1 " " i 





10 



-4 -3 -2 

10 10 10 



10 




-20 
-40 
-60 
-80 

0.2 
0.1 


-0.1 



i i 1 1 mi 1 — i i i 1 1 1 ii 1 — i i i 1 1 hi 1 — i i i 1 1 mi 1 — i r 



^ 'A 

s / 



2 2 

A = - 1 GeV 



I I I I I I I I I L 



1 1 1 1 1 1 1 III 1 1 1 1 1 1 1 If 


i i 1 1 1 1 1 1 i i i 1 1 1 1 1 1 \^\ i 


- ^ 


A 2 = -lGevV : 

— / 

S 






T 1 1 





10 



-4 -3 -2 -1 

10 10 10 10 



Figure 9: The relative NLO radiative corrections (187) and (189) are plotted versus £ for the mod- 
ulus (up) and phase (down) of S H for A 2 = (left) and A 2 = — 1 GeV 2 (right): CS (dashed) and 
MS (dash-dotted) scheme. Thick (thin) lines refer to the "hard" ("soft") gluon parameterization, 
where the scale setting prescriptions are the same as in Fig. 7. 

and is now only about ±1% or even smaller in both schemes. At this order, the renormalization 
scale dependence arises and the modulus of the CFF can vary of up to 6% [8%] by changing the 
scale in the CS [MS] scheme. We will come back to these numbers and compare them with the 
actual NNLO corrections in the CS scheme, evaluated in Sect. 6.4. We should stress here that the 
uncertainties with respect to the factorization scale setting are 'maximized' at LO. Or, in other 
words, revealing GPDs from experimental data in this approximation means that one does not 
know at which resolution scale /i 2 this information was extracted. Was it <2 2 , Q 2 /2, or . . . ? 

6.3.2 Flavor singlet sector 

For the numerical studies in the flavor singlet sector we use the same scale prescriptions and 
normalization conditions as in the preceding section. As input we alternatively take the 'soft' and 
'hard' gluonic ansatze (178) and (179), respectively. The Wilson coefficients for the CS are listed 
in Eqs. (101a) and (101b) and the evolution operator can be read off from Eq. (121). The results 
for the MS scheme are collected in Sect. 4.3. 
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Figure 10: Evolution of the flavor singlet CFFs s 7i. The moduli (up) and phases (down) are 
plotted for fixed Q 2 = 5 GeV 2 (left), Q 2 = 25 GeV 2 (right), and A 2 = -0.25 GeV 2 versus f in 
LO (dotted) and NLO for the CS (dashed) and MS (dash-dotted) scheme. The ansatze and scale 
setting prescriptions are the same as in Fig. 9. 

In Fig. 9 we plot the relative NLO corrections (187) and (189) for the moduli (up) and phase 
differences (down) at the input scale <2q = 2.5 GeV 2 , again for the two extreme values of the 
momentum transfer squared: A 2 = (left) and A 2 = — 1 GeV 2 (right). As in the nonsinglet 
case, the variation of the phase differences from LO to NLO is not large and does not exceed 
0.09tt rad for A 2 = -1 GeV 2 and is even smaller for A 2 = 0. However, the NLO corrections to 
the moduli have now a wider variety. This is related to the fact that at NLO the gluons enter 
the hard scattering part of the DVCS amplitude the first time. For the 'soft' gluon ansatz (thin 
lines) they lead to a small decrease of the CFF of about 10 — 20% and 15 — 30% for the CS and 
MS scheme, respectively. In contrast, if the gluon is 'hard' (thick lines) it cancels partly the sea 
quark dominated LO contribution and reduces so drastically the modulus of the CFF. Since in 
this scenario the gluonic part grows with decreasing £ faster than the sea quark one, the modulus 
of the CFF monotonously decreases, too. For our ansatz the reduction reaches 80% at very small 
£ and large A 2 . As in the nonsinglet case, we observe again that the radiative corrections are a 
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bit smaller at A 2 = 0. Another similarity is that in the CS scheme they are up to 5 — 20% smaller 
than in the MS one. Although the size of perturbative corrections can be very large, the phase 
differences are still small. This is caused by the fact that the phase is dominated by the leading pole 
of the ansatz. Let us stress that the rather large corrections to the hard scattering part, induced 
by gluons, should not be considered an argument against the applicability of perturbation theory. 

We study now the evolution effects to LO and NLO approximation. The ratios (193) are plotted 
in Fig. 10 versus f for two fixed values Q 2 = 5 GeV 2 and Q 2 = 25 GeV 2 , where A 2 = -0.25 GeV 2 . 
The moduli (up) in the 'soft' gluon scenario (thin lines) are relatively mildly affected by evolution 
and they grow with decreasing £. Both the modification of the LO (dotted) prediction and the 
difference between the CS (dashed) and MS (dash-dotted) schemes are rather small. However, 
we remark that the NLO contributions to the anomalous dimensions are getting large in the 
small £ region, e.g., about 100%, which is eventually compensated by the evolved NLO Wilson 
coefficients. We also found that the off-diagonal entries in the anomalous dimensions cannot be 
neglected anymore. Their contribution to the net result can grow from a few percent in the large 
£ region to over 25% in the small one. This is related to the fact that the off-diagonal entries (139) 
contain now j = poles, which arise from the LO anomalous dimensions (108) and (109). In 
the "hard" gluon scenario the evolution effects, the NLO corrections and the scheme dependence 
are quite large. The NLO corrections to the evolution are dominated by those to the anomalous 
dimensions. We again observe that the NLO corrections are smaller in the CS scheme. The 
scheme dependence partly arises from the NLO Wilson coefficients, yielding in the former scheme 
smaller corrections, which evolve with the LO evolution operator, however, also due to the off- 
diagonal part in the anomalous dimensions. Corresponding to the evolution effects that appear 
in the moduli, the phase differences in the "soft" scenario are much smaller than in the "hard" 
one. However, also in the latter case they cannot be considered large. Again, we see that at least 
within our ansatze the phase is protected from radiative corrections, since their leading pole is in 
the vicinity of j = 0. 

Table 2 lists the changes of the CFF that come from the variation of the factorization and 
renormalization scales, see Eq. (194). The first row demonstrates that the factorization scale 
variation in LO is correlated with the evolution, compare with dotted lines in Fig. (10). In the 
fixed target region the evolution and the associated variation is weak. However, approaching 
the small £ region, its strength is growing, in particular for the 'hard' gluon ansatz. One would 
normally expect that the scale variation is getting smaller at NLO. But this is not the case for the 
small £ region, rather the sign is reversed and its magnitude increases. This behavior completely 
differs from the one in the fixed target kinematics and it tells us that the factorization logs in 
the small £ region are enhanced. Hence, already from these NLO findings one might wonder 
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variation, order/ £ 


10" 5 


10" 3 


10" 1 


0.25 


0.5 


A*/, LO 


13.2 [49.5] 


9.4 [29.4] 


2.1 [7.3] 


-1.5 [1.7] 


-4.8 [-2.6] 


/i/, NLO (CS) 


-26.8[-67.6] 


-15.8 [-30.1] 


-2.8 [0.4] 


-0.5 [2.2] 


-0.3 [1.4] 


NLO (MS) 


-32.0[-48.8] 


-19.1 [-40.5] 


-3.6 [0.2] 


-0.6 [2.6] 


0.0 [2.0] 


fir, NLO (CS) 


9.3 [46.8] 


7.8 [26.8] 


6.5 [11.1] 


4.6 [7.0] 


2.4 [3.8] 


At r , NLO (MS) 


13.8 [76.5] 


12.1 [42.3] 


10.9 [17.7] 


8.7 [12.2] 


6.1 [8.0] 



Table 2: Relative changes (194) in percent of the singlet CFFs ®H within the separate variation 
of the factorization and renormalization scale from Q 2 /2- ■ -2Q 2 . Here we used the soft [hard] 
conformal GPD moments (178) [(179)] and set Q 2 = 4GeV 2 and A 2 = -0.25 GeV 2 . 

whether a perturbative treatment of the evolution in the usual manner is justified. The change of 
the renormalization scale yields variations of the few to ten percent in the 'soft' gluon scenario, 
and, with decreasing £, to much larger ones in the 'hard' gluon scenario. These numbers reflect 
simply the size of the NLO corrections for the former and latter scenario, respectively, and do 
not necessarily indicate that the perturbative expansion of the Wilson coefficients is ill-defined. 
Whether this is the case or not can be only clarified by a NNLO evaluation. 

6.4 Radiative corrections beyond NLO 

To investigate the radiative corrections in NNLO, we use the CS scheme within the Wilson coeffi- 
cients (91a)-(91c) and (101a)-(101c), as well as the evolution operators (116) and (121). We stress 
again that the perturbative expansion is consistently done as a power series in a s /2ii to the order 
N 2 LO, where the running of the coupling is described to the same order. We again calculate CFF 
H,, where the model ansatz for conformal GPD moments, scale setting, and normalization of the 
running coupling are spelled out above. We remind that the mixing term, appearing at three-loop 
level in the anomalous dimensions is unknown. Fortunately, we found that at NLO the mixing 
term in the MS scheme is small (tiny) for flavor (non-) singlet CFFs in the fixed target kinematics. 
Therefore, we expect that it is justified to neglect a NNLO mixing term for these quantities in the 
CS scheme. Unfortunately, this is not true for the singlet part at smaller values of £, where we 
observed at NLO about 30% effect at £ = 10~ 5 and Q 2 = 100 GeV 2 . Roughly speaking, we would 
presume that the mixing at NNLO is given by the contribution of the diagonal NLO anomalous 
dimensions times (— /3o)a s /2ir ~ 0.4, which is additionally suppressed by the initial condition. All 
together, we expect that for small £ the neglected mixing contributes to the net result at the 10% 
level. 
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Figure 11: The modulus (left) and the phase (right) of the rescaled singlet CFF £ ,S H(^, Q 2 ) 
versus Q 2 . Here the 'soft' gluon ansatz is used within A 2 = —0.25 GeV 2 to LO (dotted), NLO 
(dash-dotted), and NNLO (solid). 

In Fig. 11 we visualize the general features of the radiative corrections for fixed target and 
collider experiments up to NNLO in the parity even sector. Thereby, we employ the flavor singlet 
CFF ®H within the 'soft' gluon ansatz (178) at A 2 = —0.25 GeV 2 . As pointed out in the preceding 
section, in this ansatz both the modulus, scaled with £, (left panel) and the phase (right panel) 
are mildly affected by NLO corrections (dashed). As it can be seen, the NNLO corrections (solid) 
are insignificant for both of them at the input scale Q 2 = 2.5 GeV 2 over the whole £ region. As 
long as we stay away from the very small £ region, the perturbative prediction for the evolution is 
stable, starting at NLO, too. But approaching the small £ region, NNLO corrections are growing 
in size, which already shows up in a splitting of the £ = 10~ 3 NNLO and NLO trajectories, smaller 
for the modulus and larger one for the phase. The £ = 10~ 5 NNLO trajectory, compared to the 
NLO one, is affected by rather large corrections, which are of the same size as the NLO ones, 
but with a competing overall sign. With increasing Q 2 the NNLO trajectory approaches the LO 
one. Such an ill behavior reflects the competition of the Bjorken limit, i.e., Q 2 — > oo, and the 
high energy limit, i.e., l/£ — > oo. Indeed, the expansion parameter is rather m(l/£)a s (Q)/27r. 
Since the slopes of the NLO and NNLO trajectories are getting closer at large values of Q 2 , one 
can easily imagine that using a larger input scale would alleviate this problem. Indeed, it entirely 
originates from the resummed poles in the anomalous dimensions at j = and so it is universal 
(process independent) and the same one that appears for parton densities. Below we will come 
back to this issue. The changes of the phases with respect to Q 2 is for any given £ smaller than 
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Figure 12: The relative NLO (dashed) and NNLO (solid) radiative corrections to the flavor 
nonsinglet CFF NS H in the CS scheme at the input scale Q 2 = 2.5 GeV 2 (left) and Q 2 = 10 GeV 2 
(right) and A 2 = —0.25. The moduli (189) and phase differences (187) are shown in the upper 
and lower panels, respectively. The results without and with sea quarks are shown as thin and 
thick lines, respectively, which are indistinguishable for NNLO moduli. 

0.037T rad and according to Eq. (191) they are bound by the ratio of moduli. In the small £ region 
the phase approaches the value tt/2 and so the CFF is dominated by the imaginary part. This 
value is driven by both the pole, which is in our ansatz (174) in the vicinity of j = 0, and the 
essential singularity of the evolution operator at j = 0, resulting from poles in the anomalous 
dimensions (108) and (109). 

In Fig. 12 we present a more detailed view on the radiative corrections in the flavor nonsinglet 
sector for the moduli (up) and the phase differences (down). The left panels show the radiative 
corrections at the input scale Qq = 2.5 GeV 2 , while in the right panels we evolve the input to the 
scale Q 2 = 10 GeV 2 . This gives a measure of the radiative corrections arising from the evolution 
operator. Let us first discuss those for the moduli that arise at the input scale, i.e., the left 
upper panel. The NLO corrections for both a valence-like ansatz (thin) and one with a sea quark 
admixture (thick) have for the relevant fixed target kinematics a variance of about 20%. In NNLO 
this is reduced to the 5% level. More precisely, the radiative corrections reduce the CFF both 
without and with sea quark ansatz by about 2 — 5%. If we evolve them to 10 GeV 2 then the 
NNLO radiative corrections further reduce by 1 — 2%, while at NLO the variance of them is of 
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Figure 13: The relative NLO (dashed) and NNLO (solid) radiative corrections in the CS scheme 
are plotted versus £ for the modulus (up) and phase (down) of singlet 7i for A 2 = —0.25 at 
the input scale Q 2 , = 2.5 GeV 2 (left) and Q 2 = 100 GeV 2 (right). The results for the "soft" and 
"hard" gluon ansatz are shown as thin and thick lines, respectively. 



about 13%, cf. right upper panel. The radiative corrections to the phases are already dictated 
by those to the moduli. Already to NLO, they are smaller than 0.057T [0.087r] rad for the ansatz 
without [with] sea quarks and shrink further with evolution. They become tiny at NNLO, see 
lower panels. Obviously, there is an improvement in the perturbative expansion, which is roughly 
of the same order we estimated from the variation of the scales, see Tab. I. However, a closer look 
to the separate contributions, arising from different color factors, shows that there is a cancellation 
between the C\ and the /?o proportional terms [101]. The latter is negative and about two times 
larger than the former, positive one, Hence both of them partly compensate each other. Without 
this delicate cancellation the NNLO corrections would be almost on the 10% level. Let us mention 
that the scale dependencies are now almost £ independent and are of the order of —1.5% and 3% 
for the factorization and renormalization scale, respectively. 

In Fig. 13 we display the relative radiative corrections for the flavor singlet CFF in a manner 
analogous to Fig. 12. However, we evolve the quantities in questions to a scale of 100 GeV 2 , shown 
in the right panels. From the left panel, where corrections are given at the input scale, we certainly 
realize that the large negative NLO corrections (thick dashed) to the modulus in the 'hard' gluon 
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order/ £ 


10" 5 


1(T 4 


10" 3 


io- 2 


10" 1 


0.25 


0.5 


NLO 


2.4 [24.9] 


2.8 [21.0] 


3.5 [18.1] 


5.0 [15.8] 


5.8 [10.9] 


4.4 [7.1] 


2.4 [3.9] 


NNLO 


-1.6 [3.4] 


-0.6 [5.6] 


0.3 [6.5] 


0.6 [5.7] 


2.2 [6.7] 


3.5 [6.9] 


3.7 [5.9] 



Table 3: Variation (194) of ^ in percent within the change of renormalization scale from Q 2 /2 
to 2Q 2 . Here we used the soft [hard] conformal moments (178) [(179)] and set Qq = Q 2 = 4GeV 2 
and A 2 = -0.25 GeV 2 . 

scenario are shrunk to less than 10% (thick solid), in particular in the small £ region. For the 
'soft' gluon case the NNLO corrections (thin solid) are reduced to ±5%. However, for £ ~ 0.5, 
the corrections are reduced only unessentially and are still around 5% and 10% at NNLO level. 
The phase differences (lower panel) are becoming tiny at NNLO. If evolution is now switched on, 
our findings drastically change. For £ > 5 ■ 10 -2 they are stabilized for the moduli on the level 
of about 3% at Q 2 = 100 GeV 2 . However, they start to grow with decreasing £ and reach at 
£ w 10~ 5 the 20% level. It is remarkable that the relative sign of the NLO and NNLO corrections 
change and that they are becoming independent of the input. This behavior is also reflected in the 
phase differences, which decrease to — 0.027T rad. As already explained above, this breakdown of 
perturbation theory stems only from the anomalous dimensions and is thus universal, i.e., process 
independent. 

Finally, we comment on the scale dependencies. As it has been already seen in Tab. 2, the 
variation within the factorization scale increases in the small £ region with the perturbative order. 
To NNLO, we find for instance at £ = 10~ 5 a variation of 44% [105%] for the 'soft' ['hard'] 
gluon ansatz. This is about two [1.5] times larger than that observed at NLO, where the sign 
is alternating. This simply reflects the breakdown of perturbative expansion of the evolution 
operator, as we have already seen. The NLO estimates of the higher order corrections, obtained 
by the variation of the renormalization scale, were for the 'soft' gluon scenario comparable to 
the actual NNLO result at the input scale. However, the corresponding estimates for the 'hard' 
gluon ansatz in the small £ region were too pessimistic, substantially overestimating the calculated 
NNLO corrections. Since our input scale Ql = 2.5 GeV 2 in Tab. 2 was lower than the average scale 
Q 2 = 4 GeV 2 it might be that these large estimates are partly contaminated by the factorization 
logs. In Tab. 3 we show the renormalization scale dependence, but now for the input scale 
Qq = 4 GeV 2 . Compared to Tab. 2 the modifications in NLO are not large. We also realize that 
the renormalization scale independency is improved at NNLO level. 

To summarize our findings, we saw that NLO radiative corrections are moderate in the nonsin- 
glet case but can be rather large in the singlet sector for a 'hard' gluon ansatz. The factorization 
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scale dependence, substantial at LO, becomes for kinematics of fixed target experiments small 
already at NLO. However, it is getting worse in the small £ region. Interestingly, we also observe 
a scheme dependence at NLO of the order of 10% to 20%, which entirely arises due to skewness 
effects, where the radiative corrections are more pronounced in the MS scheme than in the CS 
one. At NNLO we have found that the perturbative corrections are getting reasonably small at 
the input scale and for the evolution in the fixed target kinematics. Both the factorization and 
renormalization scale dependencies are reduced to the level of a few percent. So far these findings 
suggest that the perturbative expansion is a reliable tool. However, we saw that the perturbative 
expansion of the evolution breaks down in the small £ region. Fortunately, this breakdown is 
universal, and as long as one precisely defines the scheme and the approximation, perturbation 
theory can be used as a tool to relate different processes. We will demonstrate this in the next 
section. 



7 Fitting procedure of experimental data 

The DVCS data measured in fixed target [63, 64] and collider [65, 66, 67] experiments have been 
confronted in the literature with theoretical predictions, e.g., color dipole model [168], collinear 
factorization approach within an aligned jet model inspired GPD ansatz [169] and the minimal 
dual GPD parameterization [91], see also Ref. [68] for a first analysis within the double distribution 
ansatz. Certainly, confronting model ansatze with DVCS data is a too rigid approach and should 
be considered only as a first step towards extracting GPD parameters within a given ansatz. A 
fitting procedure for DVCS data in the small £ and large Q 2 region has been proposed within 
the double log approximation of the deeply virtual Compton scattering amplitude to LO accuracy 
[166]. There a pomeron inspired GPD model was employed 20 , where the singlet quarks were 
dynamically generated. Thus, the number of fitting parameters could be reduced to three, namely, 
normalization, slope-parameter, and input scale. 

We will now demonstrate that the Mellin-Barnes representation of the CFFs is appropriate 
for a more general GPD fitting procedure. Technically, we use the standard fitting routine MI- 
NUIT [171]. This routine calls a FORTRAN code that evaluates the CFFs from conformal GPD 
moments, depending on a few fitting parameters. 

20 A Reggc pole model for the virtual Compton scattering amplitude was also used for a fitting procedure in Ref. 
[170]. 
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7.1 Setting the scene 

The DVCS amplitude interferes with the Bethe-Heitler bremsstrahlung one and so we have a rich 
selection of observables, mainly in the interference term. The decomposition in terms of CFFs is 
generally challenging. We will deal here with the easiest case, namely, the fitting of small XBj — 2£ 
data for the DVCS cross section. In this kinematics the interference term, integrated over the 
azimuthal angle, can be neglected and the DVCS cross section can be extracted from the photon 
leptoproduction one by subtracting the Bethe-Heitler bremsstrahlung cross-section. In the DVCS 
cross section we can safely neglect terms that are kinematically suppressed by £ 2 : 



2<r2 



: (W,A 2 ,Q 2 



dA 2 y ' ' ' Q 4 W 2 + Q 2 



A 2 , ,9 ~ 2 

H 



\n\ 2 — — \£\ 2 



AMI 



(e,A 2 ,Q 2 ) e Q2 . (195) 



Q 2 



2W ! + Q 

Here we have expressed the scaling variable £ in terms of the photon virtuality Q 2 and the photon- 
proton center-of-mass energy W, defined by W 2 = (Pi + qi) 2 . The leading Regge trajectory, 
appearing in Ti., arises from mesons with generic intercept a(0) ~ 1/2, which is less than the 

~ 2 

intercept «p(0) « 1 of the pomeron dominated CFFs 7i and £. Thus, the squared CFF Ti 
can be neglected, too, since it is approximately suppressed by one power of £. More care has 
to be taken about the remaining combination of H, and £ CFFs. Taking the mean value of 
<C A 2 ^>= — 0.17GeV 2 , which has been measured by the HI collaboration [67] for |A 2 | < 1 GeV 2 , 
we find that the helicity flip contribution is in the A 2 integrated cross section kinematically 
suppressed as 

<A 2 > „ _ 2 



4M 2 



10" 2 . (196) 



Hence, in this kinematical region, it might be justified to neglect the squared CFF |£| 2 . But in 
the differential cross section at larger values of —A 2 , |£| 2 might contribute to some larger extent. 
Since it is not possible to separate by the present data set the H, and £ contributions, we simplify 
our analysis by neglecting the latter one. All together, the DVCS cross section reduces in the 
small ^-region to: 

^pr.A'.fl') » ,J: a *„™* |H|' U = -$—A\tt\ ■ (197) 



dA 2 v ' ' y (2W 2 + Q 2 ) 2 W 2 +Q 2 > 1 V 2H/ 2 + Q 2 

As we have clearly pointed out in the preceding section, the perturbative expansion of the 
evolution operator is ill-defined in the small xb] — 2£ region for the (singlet) parity even sector. 
Nevertheless, we have argued that this should not affect the task of relating different processes 
within perturbative QCD. One can ask the question: How to resum this alternating series of 
ln(l/£)a s /27r terms? An answer is certainly needed for the analysis of a large amount of high 
precision data, as it is the case in DIS. Here, indeed, some progress has been recently reported, 
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see, e.g., Ref. [172]. Concerning the situation in DVCS, the solution of the problem would certainly 
improve our partonic interpretation of the nucleon content, however, is rather irrelevant for the 
analysis of present experimental measurements, as we will see. This problem will also affect the 
forward limit of conformal GPD moments, which provide the Mellin moments of parton densities. 
Since for fitting parton densities to experimental data it is also necessary to precisely specify 
the procedure, e.g., definition of perturbative expansion of evolution operator, flavor scheme, and 
running of the coupling, we do not rely on any of the standard parameterizations. Rather, we fit 
the data ourselves, within our specifications, and then compare results with a specific choice of 
parton density parameterization from literature. The perturbative expansion of the DIS structure 
function F 2 reads: 

F 2 (x Bj , Q 2 ) = —— I dj a-/ [Ql Cj (a s (Q 2 )) q 3 {Q 2 ) + Q 2 NS N %(a s (Q 2 )) N %(Q 2 )] , (198) 

where Cj = ( s Cj, Cj) and Cj are the DIS Wilson coefficients corresponding to the structure func- 
tion F%. They can be found in Ref. [141], where we set the spin label n = j + 1. Again, we equate 
the factorization and renormalization scales with the photon virtuality Q 2 and consistently com- 
bine the perturbative expansion of the Wilson coefficients with the one of the evolution operator. 
The evolution of the parton density moments is governed by the evolution equations (115) and 
(119). These moments are related to the conformal GPD ones in the forward kinematics (161), 
e.g., cf. (A.16), 



QW)=[Q H .)(ri = > A =0 >i4)- (199) 

We will employ in this kinematics our ansatz (175) and (176), see also Ref. [173], where we 

neglect for simplicity the flavor nonsinglet contribution. Moreover, instead of decomposing the 

singlet quark contributions in valence and sea quarks, we use the effective parameterization 

B(l - a s (0) + 7,8) 1 1 m . „, . . 

1 0{r] 2 ). (200) 



^ = ^ B(2 - , ; 0),8) 1 - * h .,y 



From the inspection of the standard parameterizations of parton densities, we found that the 
contributions from the flavor nonsinglet sector generally don't exceed the 10% level. Here it is 
effectively included in the parameterization of the singlet quark contribution. Pre-fitting the data 
we found that the fits are almost insensitive to the parameters a' s , a' G and that AMj; = AMq 
is the preferred value. According to the common fits to vector-meson electroproduction data in 
the small £ region, we set a'^ = a' G = 0.15 GeV -2 . Moreover, we neglect the j-dependence in the 
cut-off masses Mj 1 and M- , i.e., put AMs = AM G = 0. As relevant fitting parameters we thus 
choose 

Nv, « E (0) , Mq , N G , « G (0) , Mq . (201) 
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Obviously, in fits of the DIS structure function F 2 we have only four parameters, where «e(0) 
and «g(0) actually determine the powers of xbj and iV^ and Nq the normalization. Note that 
although in our ansatze the conformal GPD moments are normalized to H^ =1 (r], A 2 = 0, /ig) = N% 
and if ^(77, A 2 = 0, /ig) = Nq, we do not impose the momentum sum rule N-% + Nq = 1 here 21 . 
This minimal parameterization is not so well suited to provide a high quality fit to the DIS data. 
This is not our goal here; rather we would like to relate the DIS data with the DVCS data, which 
have much larger error bars. The A 2 slope of the DVCS amplitude is controlled by the cut-off 
masses and Mq. The normalization and the ^-dependence of this amplitude at A 2 = is as 
in DIS fixed by the remaining parameters N%, «s(0), Nq, «g(0). Within the assumption that the 
skewness parameter is negligible in the conformal GPD moments, we also loose the possibility to 
control the normalization of the DVCS amplitude by the skewness dependence. As demonstrated 
above by inspection of radiative corrections in different schemes, compare dashed and dash-dotted 
lines in Fig. 9, the differences in normalization, caused by the ^-dependence, is almost skewness 
independent. Within our ansatz we found a 10% — 20% effect for the modulus of the amplitude, 
which means that the cross section would differ of about 20% — 40%. The inclusion of skewness 
dependence, controlled by a corresponding parameter, will be considered somewhere else. In our 
simplified model ansatz we do not explicitly adjust the normalization of the DVCS amplitude 
relative to DIS one. However, for given mean value <C A 2 ^> the normalization of the DVCS 
amplitude, integrated over A 2 , is controlled also by the parameters and Mq , which determine 
the A 2 -slope. 

7.2 Lessons from fits 

In Fig. 14 we confront the outcome of a simultaneous x 2_ fit to the DVCS (39 data points) and 
DIS (85 data points) data in the CS scheme to NNLO accuracy. We equated the factorization and 
renormalization scale with the photon virtuality, used the conformal GPD moment ansatze at the 
input scale /ig = 4GeV 2 , fixed the number of flavors to n/ = 4, and used for a s (2.5 GeV 2 )/n = 
0.0976 22 . In the upper left panel we display the fit to the HI [67] and ZEUS [66] DVCS data versus 
—A 2 for fixed Q 2 and W, while the upper right and lower left ones show the DVCS cross section, 
integrated over |A 2 | < 1 GeV 2 , versus Q 2 and W dependence, where the remaining variable W 
or Q 2 is fixed. The fit to the DIS HI data [174] is plotted in the lower right panel, where for 

21 Note that the whole momentum fraction region, i.e., < x < 1, contributes to this rule, however, our fits 

constrain only the small XBj region, where xbj < x < 1. For technical reasons, however, we assume a certain large 

x or j behavior. A modification of this behavior, which is only weakly constrained in our fits, might be used to 

restore the momentum sum rule. An improved treatment will be given somewhere else. 

22 Employing the standard procedure for the NNLO running of coupling, we find that this value corresponds to 

a s (M|) = 0.114, at the standard reference scale M z = 91.18 GeV. 
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Figure 14: Simultaneous fit to the DVCS and DIS data in the CS scheme to NNLO. Upper 
left panel DVCS cross section for Q 2 = 4 GeV 2 and W = 71 GeV (circles, dashed) as well as 
Q 2 = 8 GeV 2 and W = 82 GeV (squares, solid) [67]. Upper right panel DVCS cross section (|A 2 | < 
1 GeV 2 ) versus Q 2 for W = 82 GeV (HI, circles, dashed) and W = 89 GeV (ZEUS, triangles, dash- 
dotted) [66]. Lower left panel DVCS cross section versus W Q 2 = 4 GeV 2 (HI, circles, dashed), 
Q 2 = 8 GeV 2 (HI, squares, solid), and Q 2 = 9.6 GeV 2 (ZEUS, triangles, dash-dotted). Lower 
right panel shows F 2 (x B j, Q 2 ) versus Q 2 for a; Bj = {8 ■ 10~ 3 , 3.2 • 1(T 3 , 1.3 • 10" 3 , 5 • 10" 4 } [174]. 
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0.116 


0.167 


1.14 


1.34 


0.535 


1.09 


1.59 


95 


0.80 


2.2 


NNLO (CS) 


0.114 


0.167 


1.14 


1.17 


0.571 


1.07 


1.39 


91 


0.77 


2.2 



Table 4: Parameters extracted from a simultaneous fit to the DVCS cross section and DIS 
structure function F 2 , where a s (M z ), a' s = a' G = 0.15 GeV -2 , and AM S = AM G = are fixed 
and f4 = 4 GeV 2 . 



clarity not all points are displayed 23 . As one realizes by eye inspection, the normalization, scale- 
and A 2 -dependency are separately well described. 

From Table 4 one can see that the quality of these simultaneous fits for the 124 data points is 
satisfying. For instance, for the NLO fit in CS scheme, where now a s (2.5GeV 2 )/7r = 0.1036, we 
get x 2 = 95, i.e., x 2 /d.o.f. = 0.8. Taking into account the NNLO order corrections, the quality 
of the fit slightly improves, i.e., x 2 /d.o.f. = 0.77. As one can realize by comparison of the last 
two rows in Tab. 4, the resulting parameters remain stable. The largest modification appears in 
the cut-off masses, which reduce from NLO to NNLO by about 13%. Thereby, contribution to x 2 
coming from the eight data points of the differential cross section (upper left panel in Fig. 14), 
denoted x^ = 2-2 (last column in Tab. 4), indicates an equally good fit to the A 2 -slope in NLO 
and NNLO. 

To understand better the double role of the cut-off masses in the fitting procedure, let us 
compare the LO parameters with the NLO ones, where their changes are much more pronounced. 
Although in LO the fit to all data yields seemingly reasonable x 2 /d.o.f. = 0.85, the A 2 -slope (de- 
termined by the small squared cut-off masses M| ~ Ml ~ 0.25 GeV 2 ) turns out to be ~ 30 GeV" 2 , 
which is incompatible with ~ 6 GeV -2 indicated by the data [67]. Or, in other words, to leading 
order accuracy = 38.5 is unacceptably large. Therefore, we conclude that the relative nor- 
malization between the DIS structure function and the integrated DVCS cross section is correctly 
reproduced in the fits by reducing the latter one by forcing the steeper A 2 -dependence. This 
also implies that at A 2 nin « the normalization of the differential DVCS cross section and the 
structure function F 2 cannot be simultaneously described within our ansatz. A separate fit to the 
DVCS data yields a slope that is getting compatible with the measured one; however, the overall 
normalization is now deteriorated and so the quality of the fit is worse, namely, x 2 /d.o.f. = 4.8. 
Certainly, this is related to the fact that we have no control over the skewness dependence of the 

23 Not shown, but used in fits, are HI F 2 data for x Bj = {1.3 ■ 1CT 2 , 5 ■ 10" 3 , 2 • 10" 3 , 8 • 10" 4 , 3.2 ■ 1CT 4 , 2 • lCT 4 , 1.3 ■ 
1CT 4 ,8- 1CT 5 }. 
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conformal GPD moments . Therefore, non-trivial skewness dependence should be introduced in 
such a way to make and consequently the normalization of the DVCS cross section, smaller. 
We remind, however, that the inclusion of the skewness dependence within the spectral repre- 
sentation of GPDs is usually done in such a way that the skewness effect leads to an increase of 
\H\. 

Let us finally compare the NLO fits in the CS and MS schemes. The fit in the latter scheme 
is compared to the former one a bit off. We find x 2 /d.o.f. = 0.92, while in the CS scheme 
X 2 /d.o.f. = 0.8. Also here the largest changes appear in the cut-off masses, in particular the gluon 
one. However, the quality of the A 2 -fit to the differential cross section is acceptable, namely, 
X 2 A 2 = 4.2, compared to %a 2 = 2-2 in the CS scheme. As discussed in Sect. 6.3 above, these 
schemes differ only in the skewness dependence of the conformal GPD moments. Instead of 
adjusting the normalization by changing the skewness dependence, it is done within our ansatz by 
increasing the slope. Hence, we must conclude that a more precise extraction of GPD parameters 
requires the inclusion of the skewness effect for A 2 = 0. 

7.3 Comparison and partonic interpretation of the results 

We would like now to confront our findings with the parton densities, as obtained from global DIS 
fits. In Fig. 15, we plot the singlet quark and gluon distributions 

(S \ i rc+ioo / E f7 \ 

4) = ^ u * (<*;) c = »• a2 = °^ < 2o2 » 

respectively, to LO (dotted), NLO for the MS (dash-dotted) and CS (dashed) as well as to NNLO 
(solid) at the input scale /i 2 , = 4GeV 2 . Note that the difference between the two schemes arises 
purely due to the skewness dependence, which affects only the DVCS fit. As argued above, we do 
not expect the outcome of, e.g., our simultaneous DVCS and DIS NLO fit, to coincide with any of 
the standard parameterizations for parton densities. Indeed, we explicitly saw this for the parton 
density fits of the HI and ZEUS collaborations. However, it turns out that the flavor singlet 
quark parton density agrees very well with the parameterization of Alekhin [175], which is plotted 
with error bands. Here the difference between the both schemes and NNLO order corrections are 
rather small. As is well known, the gluon distribution is much less constrained by a pure DIS fit. 
Here our central values lie outside the error band of the Alekhin parameterization. However, our 
error band of the NLO fit in the CS scheme would overlap with Alekhin's one. Note that we have 

24 As a side effect, arising from this ill-defined fitting task, we observed that the central value of the resulting 
parameters, in particular of the cut-off masses, are getting sensitive to the accuracy of numerics and they can vary 
inside the error bands. In NLO and NNLO fits, by increasing the numerical accuracy, we observed a variation of 
the central values only on the per mil level. 
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Figure 15: The parton densities are shown for quarks (left) and gluons (right) at the input scale 
= 4GeV 2 . The meaning of the lines are the same as before: LO (dotted), NLO for MS (dash- 
dotted) and CS (dashed), as well as NNLO (solid), the last two being indistinguishable in the 
quark case. The bands show Alekhin's NLO parameterization with errors [175]. 

used the same settings for a s as Alekhin, however, the procedures for the fits is slightly different. 
As said above, we have neglected the flavor nonsinglet contribution, took a fixed nj — 4 scheme 
for both the evolution of the running coupling and the conformal GPD moments, and assumed a 
generic j-dependence of the conformal GPD moments at the input scale. However, as we realize 
a posteriori, these simplifications are justified within the error bands. 

Although GPDs are amplitudes, it was shown that for rj = they have a probabilistic interpre- 
tation in the infinite momentum frame within the parton picture [53, 54, 55, 56]. More precisely, 
their Fourier transform with respect to the transversal degrees of freedom 

d 2 A 



H(x,b) 



-ib-A 



H(x,r) = 0, A 2 = -A 2 



(203) 



can be interpreted as parton densities that besides the longitudinal momentum fraction depend 
also on the impact parameter b. The averaged squared distance of a parton from the center of the 
nucleon might be expressed by the slope of the corresponding GPD 

Jdbb 2 H{x, b, Q 2 ^ 



(b 2 )(x, Q 2 

where the slope is defined as 
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]nH{x,ri = 0, A 2 , Q 2 
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In Fig. 16 we display the resulting slope for the flavor singlet quark combination (left) and gluon 
(right) GPD, respectively. Here we exclude the LO result, which within our ansatz fails to describe 
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Figure 16: The GPD slope, defined in Eq. (205), is shown for quarks (left) and gluons (right) at 
the input scale Q 2 = 4GeV 2 . The meaning of the lines are the same as in Fig. 15 and the band 
shows the errors for the NLO fit in the CS scheme. 



the A 2 dependence of the differentiated DVCS cross section. As observed before [166], within the 
pomeron inspired ansatz, the slope parameter slightly increases with decreasing x and is smaller 
for gluons than for quarks. As explained above, the difference between the CS (dashed) and 
MS (dash-dotted) scheme is in fact induced by our ansatz being rigid with respect to skewness. 
We would also expect that a more flexible ansatz would reduce the differences between NLO 
(dashed) and NNLO (solid) results. Compared to the LO analysis of Ref. [166], we find, e.g., for 
x = 10~ 3 and Q 2 = 4 GeV 2 , that the central value of the slope to NLO accuracy in the CS scheme 
is for quarks (gluons) about 35% (55%) smaller. In particular, for gluons is our NLO analysis 
now compatible with the slope extracted from J/ip photo- or electroproduction, see, e.g., Refs. 
[176, 177, 30]. This process is dominated by the two gluon exchange and its measured A 2 slope of 
the differential cross section, which is nearly Q 2 independent, is about 5 GeV -2 . This is consistent 
with two times the value we extracted from our DVCS analysis in the CS scheme, see right panel 
in Fig. 16. This slope yields the spatial averaged size squared for gluons of 



(£ 2 )gluon(£ 



10" 3 , Q 2 = 4 GeV 2 



0.3018;8I fm 2 (0.33^ fm 2 ) 



(206) 



to NLO (NNLO) accuracy. The central value is about 20% (10%) smaller than the one of Refs. 
[178, 179], however, still compatible within errors. We consider this a further wink that in the 
GPD phenomenology perturbative corrections should be taken into account. 
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8 Conclusions 



In this article we have derived the leading twist-two Compton form factors represented as Mellin- 
Barnes integrals in terms of conformal GPD moments with their evolution included. Thereby, we 
have used the standard framework, known from DIS, which is based on the local operator product 
expansion and dispersion relation techniques and we have confirmed known results, derived using 
other methods. For instance, this representation can be obtained in a straightforward way from the 
momentum fraction representation, which also allows other GPD related processes, e.g., the hard 
electroproduction of mesons, to be represented by a Mellin-Barnes integral. Although the original 
motivation for this representation, the solution of the LO evolution equation, is tied to conformal 
symmetry, we explicitly showed here that there is no problem in using this Mellin-Barnes integral 
representation beyond LO within the standard MS scheme, in which conformal symmetry is not 
explicitly manifested. This opens a new road for the 'global' analysis of experimental data within 
the perturbative GPD formalism to NLO accuracy. 

We have combined the Mellin-Barnes representation with conformal symmetry predictions, 
which would hold true if there existed a nontrivial fixed-point in QCD. Formallyin the perturbative 
sector one sets the (3 function to zero. Unfortunately, the return to the real (3 ^ world is plagued 
by an ambiguity, which we have shifted to the evolution equation where, consequently, it induces 
at NNLO a mixing of conformal GPD moments. We have argued that this mixing can be safely 
neglected in the fixed target kinematics, but will, however, influence the radiative corrections for 
small values of Bjorken-like scaling parameter £. 

The outcome of our numerical analysis can be summarized as follows. To leading-order ac- 
curacy, the scale setting prescription is most problematic, and this ambiguity directly translates 
into the ignorance of the scale that enters the (moments of) GPD. For fixed target kinematics, 
this problem diminishes already in NLO. There radiative corrections are moderate in the flavor 
nonsinglet sector, whereas can be larger in the singlet one, due to the appearance of gluons. With 
increasing experimental precision, one might also employ evolution to constrain the conformal 
GPD ansatze. To test the reliability of perturbation theory, we have studied NNLO corrections 
and found that they are indeed small, 5% or even less. 

In contrast to fixed target kinematics, in the small xbj region our studies have clarified the 
situation with perturbative expansion of the evolution operator being ill-defined in the kinematics 
of interest, whereas the perturbative expansion of Wilson coefficients presents no difficulties. As in 
DIS, this bad behavior arises from ln(l/£)a s terms, induced by the j = poles of the anomalous 
dimensions in the parity even sector, where j+2 is the conformal spin. Nevertheless, these poles are 
universal and, as we have demonstrated, the large fluctuation of the scaling prediction within the 
considered order does not influence the quality of fits, and, in particular, the possibility of relating 
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DVCS and DIS data. Hence, the problem of treatment or resummation of these large corrections 
is relevant primarily to our partonic interpretation of the nucleon content. As long as we precisely 
define the treatment of the evolution operator, perturbative QCD can be employed as a tool for 
analyzing data in the small Xbj region. According to this, the reported large-scale dependence 
in the hard vector-meson electroproduction might not necessarily lead to the conclusion that 
this process cannot be analyzed within the GPD formalism. We also note that progress in the 
resummation of ln(l/xBj)a s corrections in DIS has recently been achieved by two groups, see, e.g., 
Ref. [172]. 

Although we have studied only the parity even sector, one can imagine what happens in the 
parity odd one. The analytic properties of Wilson coefficients in both sectors are similar and the 
parity odd anomalous dimensions do not suffer from poles at j = 0. Hence we expect that the 
radiative corrections have similar features in the fixed target kinematics as in the parity even case 
and do not grow large in the small x-Bj region, which anyway can hardly be accessed in experiments. 

So far in the literature experimental data have been analyzed by comparing one or another 
ansatz with measurements. In particular, in the small region, where the Bjorken and the 
high-energy limit are competitors, the success of such an approach is based on trial and error. 
Indeed, as in the unpolarized DIS case, data can only be successfully fitted by a kind of fine tuning 
procedure of the relevant parameters within a clearly defined, however, ambiguous scheme. The 
most important advantage of the Mellin-Barnes representation, as demonstrated in this paper 
up to NNLO, is that it is adequate for building the fitting procedure of hard photon and meson 
electroproduction data. Certainly, a great deal of work, e.g., resummation of SO (3) partial waves, 
study of skewness dependence, and choosing an appropriate set of parameters, must be done to 
release the full power of this framework. 

To conclude, we would suggest the analysis of experimental data to NLO accuracy within the 
standard MS scheme. This drastically reduces the theoretical uncertainties, present in the LO 
approximation, particularly in relating different processes, e.g., hard photon and meson electro- 
production. The Mellin-Barnes representation allows us to write flexible fitting routines, which 
are fast enough and numerically stable, and thus it provides a reliable and systematic procedure 
to for analyzing experimental data. This we consider the main step towards a global analysis of 
experimental data, related to GPDs. 

Note added 

After the manuscript was submitted to hep-ph, the evaluation of the subtraction constant has 
been also given in Eq. (22) of Ref. [180]. Our formula (47), approximated to LO, expressed 
in momentum fraction space, and written down for CFF 7i, yields the aforementioned result 
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and contains a prescription for treating the appearing divergencies. Related work has been also 
presented in Refs. [181, 182]. 



Acknowledgements 

This project has been supported by the German Research Foundation (DFG), Croatian Ministry 
of Science, Education and Sport under the contracts no. 119-0982930-1016 and 098-0982930- 
2864, U.S. National Science Foundation under grant no. PHY-0456520, and EU project Joint 
Research Activity 5: GPDs. The authors would like to thank the theory group at the Uni- 
versity of Regensburg for its warm hospitality. D.M. is grateful for invitations at the Thomas 
Jefferson National Accelerator Facility and Service de Physique Nucleaire (Saclay). For both clar- 
ifying and inspiring discussions we are indebted to I. Anikin, H. Avakian, M. Diehl, M. Gargon, 
V. Guzey, D. Ivanov, M. Kirch, A. Manashov, B. Pire, A. Schafer, L. Schoeffel, P. Schweigert, 
L. Szymanowski, O.V. Teryaev, S. Wallon, and C. Weiss. 

A Normalization of Wilson coefficients and anomalous di- 
mensions 

Here we establish the normalization of the Wilson coefficients of the OPE (28) in the forward 
kinematics 

°C){{> = 0,Q 2 /» 2 = l,a s ) = a c)(a s ) , (A.l) 

where a cJ(a s ) and a c^(a s ) are the Wilson coefficients that appear in the perturbative description 
of the DIS unpolarized and polarized structure function F\ and g±, respectively, for fi = Q. 

The hadronic DIS tensor is related to the Compton scattering tensor (1) by the optical theorem 

= ^m7V(g, P = 2p, A = 0) , (A.2) 

where p = P\ = P2- We recall that the DIS hadronic tensor for a longitudinally polarized target 
might be written in terms of the structure functions Fi(xbj,Q 2 ), Fl{xb],Q 2 ), and gi(xBj, Q 2 )'- 

W„ v = -gl^XB^Q 2 ) + P -^F L {xB h Q 2 ) - tAe, uqp -^- gi (xBi,Q 2 ) > (A.3) 

where xbj = £ = Q 2 /2vM and A = 1 is the polarization of the target with respect to the p 
direction. Both Fi and g\ structure functions are expressed to leading power accuracy as 

CJ^-U^^W&S-^W^ 1 (A ' 4) 
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where the flavor sum runs over quarks, antiquarks, and gluons. Here (A)c a are the partonic cross 
sections, depending on the factorization scale fi. The polarized and unpolarized (anti-)parton 
distributions are denoted as q a and Aq a [(A)q a = (A)g„], respectively, and Q a are the fractional 
electrical charges 

Qu = Qc = -z , Qd = Q s = -z , Qg = — & for g luons ' ( A - 5 ) 

O O TX f 

J a=u,d,--- 

while the coefficient functions are normalized as follows 

(A) a c = 5(1 - x) + 0{a s ) for a = {u, % d, ■ ■ ■ } and (A) G c = 0(a s ) . (A.6) 

Since the evolution will yield a mixing of the quark singlet and gluon parton densities, it 
is appropriate to introduce a group theoretical decomposition in flavor nonsinglet and singlet 
densities. For three active light quarks these combinations are 

(A) NS g = 2{A) u q + 2(A)^ - (A) d q - (A)^ - (A) s q - (A)^ (A.7) 
(Afq = (A) u q + (A)^ + ( A) d q + ( A)^ + ( A) s q + (A)% , 

while for four active quarks they read: 

(A) NS g = {A)\ + (A)^ - (A) d q - (A)1? - (A) s q - (A)^ + (A) c q + (A) c q (A.8) 
(Afq = (A) u q + (A)^ + {A) d q + (A)^ + {A) s q + (A)^ + (A) c q + {Afq , 

The corresponding squared charge factors are defined in Eq. (27). In the singlet case the parton 
densities with a = {+, — } are eigenvectors of the evolution equation, i.e., all distributions satisfy 
the evolution equation 

^A-(Arq(x,^)= f^-{A) a p(-^ s ^)]{Arq{y^ 2 ) for a = {NS,+,-} (A.9) 

a/r J x y \y J 

where (A) a P are the splitting kernels. The rotation to this basis results into the replacement of 
the flavor sum in (A. 4). 

E =* E • (a.io) 

a=u,u,-- ,G a=NS,+,- 

Using the optical theorem (A. 2) and the parameterization of the Compton tensor (3) and the 
DIS hadronic tensor (A. 3) we find that the structure functions are related to the CFFs by 

C I }(^ = s 9m {s} ({= ^- i0 -- ) L- (A - U) 
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Forming Mellin moments, we find by means of Eq. (36) 



f dxx>{ Fl X(x,Q 2 )= 1 - Ql \Ia\mQ)) {If} (r/ = 0,A 2 = 0,Q 2 ),(A.12) 

^° I 9l ) a=NS,± ^ ' 3 ^ ' 3 

and, on the other hand, in terms the parton densities and coefficient moments of Eq. (A. 4) 

fdxxi\ Fl ){ X ^) = \ £ Ql («.(<?)) {2a a ) (Q 2 )> (A.13) 



a=NS,± ^ J 

with definitions 

{A)\{Q 2 ) = f dxxt(A) a q(x,Q 2 ), {A) a Cj (a s ) = [ dxx j (A) a c(x, 1, a.) . (A.14) 
Jo Jo 

From this we obtain the equalities 

a Cj(a s ) \{Q 2 ) = a cJ(a s ) a H^ = 0, A 2 = 0, Q 2 ) , (A.15) 
A a Cj (a s ) A^(Q 2 ) = a Cj A (a s ) ^(77 = 0, A 2 = 0,Q 2 ) . 

The normalization of the conformal operators (21) and (23) is chosen so that their reduced matrix 
elements in the forward kinematics reduce to those used in DIS. Taking into account that Eqs. 
(31) and (32) relate the forward matrix elements to the conformal moments Hj and Hj, we find, 
for instance, in the unpolarized case for odd j 

a=o = ^+T<Pl^Jb> = ^+T<Pl^(°)7+ (* D + ) 3 Va(0)b) = \ + % for a = {u, d, s} , 

lrM"V]\p) = -2 FT <p|^(0)A 8 7+ (< tf(0)|p> = % for a = {NS, £} , 



A=0 P| +1 " 1 J " ' 2p 



-1 (Pl^b) = TTh (pl^/(0) (< 5+Y~ G, + (0)b) = % . (A.16) 
Hence, from the identity (A.15) we find the desired normalization for the Wilson coefficients 

a cJ(a s ) = a Cj (a s ) and a cf(a s )=A a Cj (a s ). 

Taking the forward limit in the evolution equation (81) and comparing it with the moments of 
Eq. (A. 9), we also establish from the identities (A.16) the definition of the anomalous dimensions 
in terms of moments of the splitting kernel 

.1 1 r l 1 

dxx^ a P(x,a s ) = -- a ^(a s ), / dxx? A a P(x,a s ) = -- a 1 fM. (A.17) 
2 Jo 2 

Note that within our definitions, the momentum sum rule for j = 1 is established by 

±(P 2 , ftl^Pi, ft) + -^(P 2 , SjI^IA, ft) = ^ <P 2 , ft|e ++ |P!, ft) , (A.18) 
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which results in the scale independent expectation value of the energy-momentum tensor 

©++ = ^7+ D + ^ + GfG; + , (A. 19) 

projected on the plus light-cone components. In the forward case the sum rule (A. 18) reduces to 

VaO + G <iM = ^^(j>\e ++ \p) = i. (A.20) 

As a consequence of the scale independence the anomalous dimensions satisfy the relation: 

s V(«.) + G V(«.) = GG 7r(«.) + EG 7r(«.) = o. (A.2i) 



B Bases in the flavor singlet sector 

The transformation of the gluon and quark singlet conformal operators to the basis of the ± ones 
is defined in Eq. (24). In the following we drop the superscript I £ {V, A}. By making use of the 
evolution equation it can be easily shown that the two dimensional anomalous dimension matrix 
is diagonalized by the rotation 

h U \ _ TT h h ITT \-l _ / „ U TT \ (TT \~ l 



-, bJ ^{^ ^)v>r-^v l){ u,r. ( B,) 

Let us suppose that for a given scale /x the inhomogeneous term vanishes, i.e., ^j^U\^ =fl0 = 0. 
At this reference point the eigenvalues of the anomalous dimension matrix are 

% = ~ (H + GG 7, ± V / ( E ^- GG 7,) 2 + 4SG T, G ^) > (B.2) 
and the rotation matrix reads 



GG 7,'-~7, 



£/>o,/^o) = ( ^ ^ ^ |(a,W). (B.3) 

The diagonalization of the evolution equation at an arbitrary scale can be now obtained by the 
use of evolution 

Uj(ji, fio) = Ej(ji, fj^Ujijio, no)£~ l {ii, ii ) (B.4) 
where £ _1 is the inverse of the evolution operator (120) and 

Ej(fi,fi ) = \ J r , _ i | (B.5) 

v ° ex p^-/; ^(« s (^)) 

is the evolution operator in the {+, — } basis. 
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C Momentum fraction representation versus conformal mo- 
ments 



In the momentum fraction representation the Compton form factors are represented as convolution 
of the coefficient function with the corresponding GPD. In the singlet sector in which quark i^O) 
and gluon ( g O) operators mix under renormalization, we might introduce the vector notation: 

j C(x/^Q 2 /^,a s (ri\OF(x, V = ^A 2 ^ 2 ). (C.l) 

Here the column vector 

f =(g£)» F = {H,E,H,E} (C.2) 

contains the GPDs, and the row one, defined as C = \C, (1/^)^), consists of the hard scattering 
part that to LO accuracy reads 

^C(x/£,Q 2 /fi 2 ,a s (ri\0= L_l_ ie ^)+0(a s ). (C.3) 

We remark that the £ dependence in and ^ enters only via the ratio x/^. Note also that the 
w-channel contribution in the quark entry (C.3) has been reabsorbed into the symmetrized quark 
singlet distribution 

E F(*,77,A 2 ,/, 2 ) = [ q F{x^A 2 ^ 2 )T q F{-x^A 2 ^ 2 )] . (C.4) 

q=u,d,--- 

Here the second term in the square brackets with — (+)-sign for H, E (H, E)-type GPDs is for 
x > 7] related to the s-channel exchange of an antiquark. The gluon GPDs have definite symmetry 
property under the exchange of x — > —x: G H and E are even, while G H and G E are odd. 

C.l Evaluation of conformal moments 

The convolution formula (C.I) has already at LO the disadvantage that it contains a singularity 
at the cross-over point between the central region (—77 < x < 77) and the outer region (77 < x < 1), 
i.e., for x = £ = 77. Its treatment is defined by the ie prescription, coming from the Feynman 
propagator. The GPD is considered smooth at this point, but will generally not be holomorphic 
[88]. The fact that both regions are dual to each other, up to a so-called D-term contribution 
[120], makes the numerical treatment even more complicated. This motivated our development of 
a more suitable formalism in [98]. 
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To make contact with the conformal OPE, we expand the hard-scattering amplitude in terms 
of Gegenbauer polynomials with indices 3/2 and 5/2 for quarks and gluons, respectively, and 
introduce the conformal GPD moments, which formally leads to 

oo 

A 2 , Q 2 ) = 2 r^C^QV/i 2 , aM) Ftf, A 2 , / u 2 ). (C.5) 

3=0 

The expansion coefficients Cj can be calculated by the projection: 



,2/. 2 „ 2 i+1 r(i + V2) 



r(3/2)r( J + 4) 



x-/ dxCCx,^, 05.(^)1^=1) w " n 13 r 212 < /2 (x). (C.6) 
^J-i \ 3[l-x z J z C j l 1 

Note that we have here rescaled the integration variable with respect to £ and that the integral 
runs only over the rescaled central region. The conformal moments of the singlet GPDs are defined 

as 



2^r(j + 3/2) 2 7_, I o (!/)1cS ill 1 



Here j is an odd (even) non-negative integer for the (axial-)vector case. 

In order to make use of the NLO MS results given in the momentum fraction representation, 
for quark part we determine the corresponding conformal moments using the results from Ref. 
[137]. In particular, the conformal moments (128) for the axial- vector case can be read off from 
Eq. (3.32) in Ref. [137], where the normalization factor (2j + 3)/(j + + 2) must be removed 
and the remaining expression multiplied by the color factor Cf- The conformal moments (127) for 
the vector case one can easily recover from the hard-scattering amplitude given in Ref. [154], Eqs. 
(14) and (15). Thereby one has only to evaluate the difference between the vector and axial-vector 
results 

E C V(1) (x, Q/V 2 ) - Q/H 2 ) = -— ln(l - x) • (C.8) 

x 

Employing Tab. 8 in App. C of Ref. [137], we find, after removing the normalization factor 
2(2j + 3)/(j + l)(j + 2), that the difference for the conformal moments is Cf/(J + l)(j + 2). 



C.2 Conformal moments for the gluon part 

Following the method for computing moments with respect to conformal partial waves (with the 
index k) in the quark sector that was explained in detail in App. C of Ref. [137], we here give the 
necessary results for the gluon sector. 
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We introduce the notation similar to the one used in Ref. [137] 

(G(x))<£!*> = f dx G(x) X2{ l~ /2 X? Cl\(2x - 1) , (C.9) 
Jo N k-i 

with 

5/2 _ k(k + 3) 3/2 _ (fc + l)(fc + 2) 

~ ^Vp~ * ~ m + 3) ' ( } 

It follows trivially that (G(l — x)}flf = (— and one can easily make the corre- 
spondence to definitions of conformal moments given in (C.6): 

G 2 fc+1 r(fc + 5/2) 5/2 /Gr yo^ iN\(6/2) (Tin 

Cfc - r(3/2)r(fc + 4) 48 N ^ { c{2x - 1))k ~i (C11) 



It is convenient to use the following expression for the Gegenbauer polynomials: 

x l+2 . (C.12) 



x 2 (l~x) 2 r s/2 (2r n _ f nfc -i 12(2fc + 3) ^ , fk + l\ (k + i + 1 
The evaluation of the conformal moments, i.e., in our case the evaluation of the expressions 



x I fc _ 1 k(k 



i=0 

and 



o 



i=0 



consists then in calculating the Mellin moments and performing the summation. The Mellin 
moments for the functions we encounter and most of the nontrivial sums we are left with can be 
found in [150]. 

In Table 5 we summarize the conformal moments of the functions relevant to our calculation. 
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